Re normalization- group study of Anderson and Kondo impurities in gapless Fermi 

systems 
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Thermodynamic properties are presented for four magnetic impurity models describing delocalized 
fermions scattering from a localized orbital at an energy-dependent rate T(e) which vanishes pre- 
cisely at the Fermi level, e = 0. Specifically, it is assumed that for small |e|, T(e) oc |e| r with r > 0. 
The cases r = 1 and r = 2 describe dilute magnetic impurities in unconventional (d- and p-wave) 
superconductors, "flux phases" of the two-dimensional electron gas, and certain zero-gap semicon- 
ductors. For the nondegenerate Anderson model, the main effects of the depression of the low-energy 
scattering rate are the suppression of mixed valence in favor of local-moment behavior, and a marked 
reduction in the exchange coupling on entry to the local-moment regime, with a consequent narrow- 
ing of the range of parameters within which the impurity spin becomes Kondo-screened. The precise 
relationship between the Anderson model and the exactly screened Kondo model with power-law 
exchange is examined. The intermediate-coupling fixed point identified in the latter model by With- 
off and Fradkin (WF) is shown to have clear signatures both in the thermodynamic properties and 
in the local magnetic response of the impurity. The underscreened, impurity-spin-one Kondo model 
and the overscreened, two-channel Kondo model both exhibit a conditionally stable intermediate- 
coupling fixed point in addition to unstable fixed points of the WF type. In all four models, the 
presence or absence of particle-hole symmetry plays a crucial role in determining the physics both 
at strong coupling and in the vicinity of the WF transition. These results are obtained using an 
extension of Wilson's numerical renormalization-group technique to treat energy-dependent scatter- 
ing. The strong- and weak-coupling fixed points of each model are identified and their stability is 
analyzed. Algebraic expressions are derived for the fixed-point thermodynamic properties, and for 
low-temperature corrections about the stable fixed points. Numerical data are presented confirming 
the algebraic results, identifying and characterizing intermediate-coupling (non-Fermi-liquid) fixed 
points, and exploring temperature-driven crossovers between different physical regimes. 



PACS numbers: 72.15.Qm, 75.20.Hr 
I. INTRODUCTION 

In conventional metallic systems, it is well understood 
how many-body correlations induced by dilute magnetic 
impurities in an otherwise noninteracting conduction 
band can at low temperatures effectively quench all spin 
fluctuations on each impurity siteu This, the Kondo ef- 
fect, depends critically on the presence of fermionic ex- 
citations down to arbitrarily small energy scales. The 
impurity properties are sensitive to the density of elec- 
tronic states p(e) only through its value at the Fermi 
level, e = 0. Other details of the band shape have negli- 
gible effect on thedossL-temperature physics. 

A growing bodyotlll of theoretical work shows that the 
standard picture of the Kondo effect must be fundamen- 
tally revised in order to treat "gapless" systems, in which 
the effective density of states vanishes precisely at e = 
but is nonzero everywhere else in the vicinity of the Fermi 
energy. The goal of the present paper is to extend the 
understanding of this issue through a comprehensive ac- 
count of the different physical regimes exhibited by mag- 
netic impurities in gapless host materials, including de- 
tailed calculations of thermodynamic properties. 

Gaplessness may be realized in a number of physical 
systems: (1) The quasiparticle density of states in an 



unconventional superconductor can vaxy like |e| or |e| 2 
near line or point nodes in the gap.EiJ Heavy-fermion 
and cuprate superconductors are strong candidates for 
this behavior. (2) The valence and conduction bands 
of certain semiconductors touch in such a way that, for 
small |e|, p(e) is proportional to |e| rf_1 in d spatial dimcp. 
sions. Examples include PbTe-SnTe heterojunctionsji^l 
and the ternary compounds Pbi-^Sn^Se, Pbi-^Sn^Te, 
and Hgi^ajCd^T^-jeach at a temperature-dependent crit- 
ical composition.t3 (Zero-gap mercury cadmium telluride 
has been proposed as the basis for a giantjpagnetoresis- 
tance read-head for high-density storage.li-3) (3) Various 
two-dimensional electron systems — including graphite 
sheets,E-3 "flux phases" in a strong magnetic field and 
exotic phases of the Hubbard modeO — are predicted 
to exhibit a linear pseudogap. It is a matter of ongoing 
debate whether this pseudogap survives the presence of 
disorder .Ea (4) The single-particle density of states in the 
one-dimensional Luttinger model varies like |e| 2a , where 
a changes cotitinuously with the strength of the bulk 
interactions .La In all these examples, the effective den- 
sity of states can be approximated near the Fermi level 
by a power law, p(e) oc |e| r with r > 0. 

The first theoretical study of magnetic impurities in 
gapless Fermi systems was carried out by Withoff and 
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Fradkinp who assumed a simplified density of states in 
which the power-law variation extends across an entire 
band of halfwidth D: 



P(e) 




\e\<D; 
otherwise. 



(1.1) 



Using poor-man's scaling for the spin-^ (i.e., impurity 
degeneracy N = 2) Kondo model and a large- N treat- 
ment of the Coqblin-Schrieffer model — both methods 
being valid for < r < i — these authors demonstrated 
that the Kondo effect takes place only if the dimension- 
less, antiferromagnetic electron-impurity exchange poJ 
exceeds a critical value, poJ c ~ r; for J < J c , the deple- 
tion of low-energy conduction states causes the impurity 
to decouple from the band at low temperatures. 

Subsequent work has analyzed values of r Wito 2 and 
beyond. Large- N methods have been appliedtffl to mod- 
els describing magnetic impurities in unconventional su- 
perconductors, in which the power-law variation of p(e) 
is restricted to a region |e| < A. These studies, which 
may be directly relevant for Ni-doping experimentsEj on 
YBa2Cu307_a, have yielded results in general agreement 
with Ref. ^. The logarithmic dependences on tempera- 
ture T and frequency to which characterize the standard 
(r = 0) Kondo effect are replaced for r > by power 
laws. In the specific case.r — 1, these power laws acquire 
logarithmic corrections .□ For r < 1 or N = 2, any finite 
impurity concentration produces a-small infilling of the 
pseudogap which drives J c to zero. til 

Numerical renormalization-group (RG) calculations' 
for the case N = 2 have revealed a number of addi- 
tional features. At particle-hole symmetry, the critical 
coupling J c is infinite for all r > g, while for r < \ the 
strong-coupling limit exhibits anomalous properties, in- 
cluding values of the impurity entropy and the effective 
impurity moment which are noftzero (even at T — 0) and 
which vary continuously with r.B Away frorajiarticle-hole 
symmetry, the picture is markedly different .0 Progressive 
introduction of band asymmetry or of impurity potential 
scattering initially drives J c for r > | back down towards 
the large- N value paJ c ~ r; eventually, though, further 
increasing the asymmetry tends to freeze the motion of 
conduction electrons near the impurity site, leading to an 
upturn in J c . For J > J c and r > 0, the impurity entropy 
and the effective impurity moment both approach zero 
at T = 0. An electron phase shift of ir suggests that-the 
impurity contribution to the resistivity also vanishes,0 in- 
stead of taking its maximum possible value as it does in 
the standard Kondo effect .El 

The spin-s Kondo model presupposes the existence 
of a local moment at the impurity site, i.e., an im- 
purity level having an average occupancy (rid) ~ 2s. 
The more-fundamental Anderson model allows for real 
charge fluctuations on the impurity site. In the nondc- 
generate (N — 2) version of this model, mixed-valence 
(0 < (rid) < 1) and empty-impurity ({rid} ~ 0) regimes 
compete with s = \ local-moment behavior. Poor-man's 



scaling has been appliecB to an Anderson impurity lying 
inside a power-law pseudogap. The reduction in the den- 
sity of states near the Fermi level has three main effects, 
each of which grows more pronounced as r increases: the 
mixed- valence region of parameter space shrinks, and for 
r > 1 disappears altogether; there is a compensating 
expansion of the local-moment regime and, to a lesser 
extent, of the empty-impurity regime; and the value of 
the Kondo exchange J on entry to the local-moment 
regime is reduced. Since the threshold exchange for the 
Kondo effect (J c defined above) rises with r, these re- 
sults imply — at least in the cases of greatest interest, 
r = 1 and 2 — that over a large region of phase space, 
the low-temperature state possesses an uncompensated 
local moment. This should be contrasted with systems 
having a regular density of states, in which an Anderson 
impurity is always quenched at zero temper ature.EI 

This paper contains a detailed study of four mod- 
els describing a magnetic impurity in a gapless host: 
the nondegenerate Anderson model, and three variants 
of the Kondo model, representing (in the_ nomencla- 
ture introduced by Nozieres and Blandintil) "exactly 
screened," "underscreened," and "overscreened" impu- 
rity spins. The nonperturbative RG formalism originally, 
developed to describe magnetic impurities in metalsc3'E3 
is extended to provide algebraic results for the stability 
of, and properties near, the weak- and strong-coupling 
fixed points of each model. Numerical implementation of 
the RG scheme enables characterization of the thermo- 
dynamic properties at intermediate-coupling (non- Fermi- 
liquid) fixed points, and allows the study of temperature- 
driven crossovers between fixed-point regimes. 

Three effects of the pseudogap are found to be com- 
mon to all the models. (1) Over a finite fraction of pa- 
rameter space, the impurity becomes asymptotically free 
in the limit T — > 0. The size of this weak-coupling re- 
gion grows with increasing r. (2) The presence or ab- 
sence of particle-hole symmetry plays a crucial role in 
determining the low-energy behavior. It turns out that 
each model can exhibit two distinct fixed points of the 
Withoff-Fradkin type: one preserving, and the other vio- 
lating, particlc-holc symmetry. (These fixed points coex- 
ist only over a limited range of r values.) The strong- 
coupling physics is also very sensitive to particle-hole 
(a)symmetry. (3) Power-law dependences of physical 
quantities on temperature and frequency are generally 
different for sublinear and superlinear densities of states. 
Specifically, in many places where the exponent r enters 
physical quantities for r < 1, it is replaced for r > 1 by 
either 1 or r — 1. The case r = 1, of particular interest 
in the contexts of high-T c superconductivity and of two- 
dimensional flux phases, exhibits logarithmic corrections 
to simple power laws. 

Extensive results are provided for the Anderson model. 
We investigate in a systematic fashion the nature of the 
phase diagram at a fixed, positive value of r, and study 
the various trends produced by increasing r. For r > 0.5, 
it proves impossible to observe Kondo screening of an 
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Anderson impurity in a system having a pure power-law 
density of states. The suppression of the Kondo effect 
becomes less dramatic, however, if the power law is re- 
stricted to a pseudogap of halfwidth A -c D. (A pre- 
liminary version of these results was used to support the 
perturbative scaling theory presented in Ref. ||) 

Recently, Bulla et alh2 have also applied the numerical 
RG approach to the Anderson model with a pure power- 
law scattering rate, limited to cases of strict particle-hole 
symmetry. With one minor exception, the weak- and 
strong-coupling thermodynamics reported for r < i are 
in agreement with Ref. ^ and the present work. The au- 
thors of Ref. |l^ interpret these results, and certain non- 
integer exponents describing their numerical data for the 
impurity spectral function, as evidence for non-Fermi- 
liquid behavior. We demonstrate, however, that the 
fixed-point properties are precisely those expected for a 
nonintcracting gapless system. Hence, we argue that the 
weak- and strong-coupling limits can be described within 
a generalized Fermi-liquid framework. 

Previous studies of the exactly screened s = i Kondo 
model with power-law exchange have focused on the ex- 
istence and position J c of the intermediate-coupling fixed 
point, and on the thermodynamics in the weak- and 
strong-coupling regimes. Here we concentrate instead on 
the properties of the J c fixed point, which is shown to 
have a clear signature in the impurity contribution to 
total thermodynamic quantities and in a static response 
function which probes the local behavior at the impu- 
rity site. We also examine in some detail the relation- 
ship between the Kondo and Anderson models in gapless 
hosts, and conclude that the models are independent to 
a greater extent than in the standard case r = 0. 

Our investigation of the underscreened, s = 1 Kondo 
model and the overscreened, s — ^, two-channel Kondo 
model focuses on the fixed-point physics. Over a range 
of exponents < r < 1/4, each problem exhibits an 
unstable fixed point of the Withoff-Fradkin type at a 
critical coupling poJ c f=a r. The novel feature, however, 
is the existence at some exchange J* > J c of a second 
intermediate-coupling fixed point which is locally stable 
with respect to perturbations in J. The J* fixed point of 
the underscreened problem has no counterpart in metals, 
but that of the overscreened model is the generalization 
to r > of the non-Fepni-liquid fixed point identified by 
Nozieres and BlandinH For r > 1/4, the J* fixed point 
disappears, and the J c fixed point can be reached (and 
hence the Kondo effect observed) only under conditions 
of strong particle-hole asymmetry. 

Before proceeding further, we remark on a matter 
of terminology. It will be shown in the next section 
that the conduction-band density of states p(e) and the 
energy-dependent hybridization t(e) (describing hopping 
between a magnetic level and the conduction band) en- 
ter the Anderson impurity problem only in combination, 
through the scattering rate 



The exchange and potential scattering in the Kondo 
model have the same energy-dependence as T(e). In gap- 
less syst ems is is natural to assume that p(e) is given by 
Eq. (1.1) while t(e) is essentially constant. However, the 
separate forms of the density of states and the hybridiza- 
tion are unimportant provided that one is interested only 
in impurity properties. In the remainder of the paper, we 
shall therefore refer to a power-law scattering rate or ex- 
change. We shall focus mainly on the simplest case, that 
of pure power-law scattering, 



r(e) = 



role/Dp 
0, 



H < D ; 

otherwise. 



(1.3) 



T(e)=7rp(e)t 2 (e). 



(1.2) 



However, we will examine the effect of including more 
realistic features such as band asymmetry and restriction 
of the power-law variation to a finite pseudogap region. 

The organization of this paper is as follows: In Sec- 
tion [n] we describe the generalization of Wilson's numer- 
ical RG method to handle magnetic impurity problems 
with an energy-dependent impurity scattering rate. The 
three sections that follow develop the analytical aspects 
of the technique in the specific context of a pure power- 
law scattering rate, as given in Eq. (1.3). Section |l| 
deals with the discretized conduction-band Hamiltoni- 
ans that lie at the heart of Wilson's method. Section |^ 
addresses the stability of the weak- and strong-coupling 
fixed points of the four magnetic impurity models of in- 
terest, while Section ^ focuses on their thermodynamic 
properties. The reader who is already familiar with the 
models we study and who wishes to pass over the tech- 
nical details of our treatment may wish to jump directly 
to Section VI, where detailed numerical result s ar e pre- 
sented. The results are summarized in Section VII. Two 
appendices contain mainly technical details. 



II. GENERALIZED FORMULATION OF THE 
NUMERICAL RG METHOD 

In this section, we describe a generalizatinu. of Wil- 
son's nonperturbative numerical RG mcthodESO to treat 
impurity models in which the scattering rate of conduc- 
tion electrons from the impurity site is energy-dependent. 
The generalization, presented here in the context of the 
single-impurity Kondo and Anderson Hamiltonians, was 
developed independently by several groups. It, has been 
applied to the two-impurity Anderson, model J2j the two- 
impurity, two-channel Kondo modeled the Anderson lat- 
tice in infinite spatial dimensip|BSjt3 and to the single- 
impurity KondcQ and Andersonl3ll3 models with a power- 
law scattering rate. Only the last of the papers cited 
reports any technical details. Here we provide a compre- 
hensive explanation of the method. 

An alternative (but closely related) generalization of 
the numerical RG method, developed by Chen and 
Jayaprakash,E3 has been used to obtain equivalent phys- 
ical results for the Kondo model with a pure power-law 
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scattering rateJa The relationship between the two for- 
mulations is discussed in Ref. fLOl 



A. The Anderson impurity model 

The nondegenerate Anderson HamiltoniariH for a sin- 
gle magnetic impurity in a nonmagnetic host can be writ- 
ten as the sum of conduction-band, impurity, and hy- 
bridization terms: 



where 



7~La =TC c + Tid + fit, 



where 



Cko 



U d = tdn-d + Un d ^ridi, 



(2.1) 

(2.2a) 

(2.2b) 
(2.2c) 



The energies £k and td of electrons in the conduction 
band and in the localized impurity state, respectively, 
are measured from the Fermi energy; Nq is the number 
of unit cells in the host; rid = n<v\ + n-di is the total occu- 
pancy of the impurity level; and U > is the Coulomb 
repulsion between a pair of localized electrons. Without 
loss of generality, the hybridization matrix elements ik 
between localized and conduction states can be taken to 
be real and non-negative. Throughout the paper, sum- 
mation over repeated spin indices (cr in the equations 
above) is implied. 

For simplicity, we consider a spatially isotropic prob- 
lem, i.e., one in which et = eiki and ik = i(e|k|)j so 
that the impurity interacts only with s-wave conduction 
states centered on the impurity site. The energies e of 
such s-wave states are assumed to be distributed over the 
range — (l+fi)D < e < (l—fj,)D. It proves convenient to 
work with a dimcnsionless energy scale, e = e/ D. Then, 
dropping the kine tic energy of all non-s-wave conduction 
states, Eqs. (2.2) can be transformed to the following 
one-dimensional form: 



(2.3a) 
(2.3b) 



de e c\ a c ea , 

-(i+m) 

H d = tdn-d + Un^ridi, 
H h = ( " de ^fp{eD)D t(eD) (c^d, + H.c). (2.3c) 

The operator c ea , which annihilates an electron in an 
s-wave state of energy e, satisfies the anticommutation 
relations {c\ a , c e v<} = 8(e — e') 6 a ,cr'- 

In this model, the impurity couples to a unique linear 
combination of s-wave conduction states associated with 
an operator 



f0a = F~ 



de w(e) c e<T , 



(2.4) 



F 2 = 



1— /U 



(i+m) 



de w 2 (e). 



(2.5) 



The weighting function w(e) entering Eqs. (2.4) and (2.5) 
is 



U! 



(e)=yr(eD)/r 0l 



(2.6) 



where T(e) is dchned in Eq. (L2) and Tq is a reference 
value of the scattering rate (for example, that at the 
Fermi level). With these definitions, the hybridization 
term in the Hamiltonian can be rewritten 

H h = v/TVDA F (fld a + H.c). (2.7) 



B. The Kondo impurity model 

The Kondo model!! describes the interaction between 
a conduction band and a localized impurity which has a 
spin s of magnitude \J s(s + 1). The Hamiltonian is 



(2.8) 



wher e H . c is the conduction-band Hamiltonian given in 



Eq. (|2.2|a) and 



k.k' 



Jk,] 



N 



-ko-2' 



Vi 



k,k' f 



N 



(2.9) 



The first and second terms in H. s describe exchange and 
potential scattering, respectively. 

The impurity-spin- i version of the Kondo model can 
be regarded as a li miti ng case of the nondegenerate An- 
derson model [Eq. (§j|)]. If-e^, U+e d > T, k B T (where 
ksT is the thermal energy scale), then single occupancy 
of the Anderson impurity level is overwhelmingly favored 
over zero or double occupancy, in effect localizing a pure- 
spin degree of freedom at the impurity site. The exchange 
and potential scattering coefficients can Jse determined 
using the Schrieffer- Wolff transformation:' 



Jk,k> 
Vk,k' 



1 



1 



Md\ \U + e d \ 
1/1 1 



2VM \U + e d \ 



ik*k' , 
£k^k' 



(2.10a) 
(2.10b) 



Equations ( 2.10 ) imply that the exchange and poten- 
tial scattering both exhibit the same dependence on k, 
and hence (in a spatially isotropic problem) on e. Thus, 
just as for the Anderson model, the Kondo impurity in- 
teracts with a si ngle linear combination of conduction 
states. Equation (2J3) can be rewritten 
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1-/1 



Hs = I de v / p(eD)D / de' v / p{s'D)D x 

[J( £ £>, e'D) ow ■ s + V(eD, e'D) 6^,] 4 a c e>!T , 
= D [p Q J Q \rT aa , + poVoS^] F 2 flf 0a ,, (2.11) 

where fo a and F are denned in Eqs. ( |2.4| )"(^6|); 
PqJo and poVo are reference values of p(e)p(e')J(e, e') 
and \] p(e)p(e')V(e, e'), respectively. 

While we shall primarily focus on the conventional 
Kondo model [Eq. ( |2.8| ) with s — |], we shall also present 
results for the s = 1 model and for the s = 4, two- 
channel model. Following Ref. ^l], we refer to these three 
variants as the "exactly screened," "underscreened," and 
"overscreened" cases, respectively. r— . 

The A^-channel Kondo Hamiltonian,Eil describing an 
impurity spin degree of freedom interacting with N c > 1 
degenerate bands (or "channels") of conduction elec- 
trons, corresponds to Eq.(|2.8[) with 



He 



E 



ekC kj"cr Ck i cr ' 



(2.12) 



and 

k,k'j L 



J, 



(J) 



k.k' t 1 

C kjcr 2 <Tcrcr ' Ck 'i '' 



v. 



k,k' f 

C kj<T Ck 'i CT 



(2.13) 



where j = 1, 2, . . . , N c is the channel index. In this paper 
we treat only the channel-symmetric version of the two- 
channel problem, i.e., we take N c — 2, J^, = J^u') an< ^ 

V kl = V ul- (Multichannel variants of the Anderson 
model also exist, but they lie beyond the scope of the 
present work.) 



C. Tridiagonalization of the conduction- band 
Hamiltonian 

Given the form of the weighting function w(e) which 
defines fg a — the particular linear combination of de- 
localized states that interacts with the impurity de- 
grees of freedom in the Anderson or Kondo model - 
the conduction-band Hamiltonian— can be mapped ex- 
actly, using the Lanczos procedure,^ onto a tight-binding 
Hamiltonian describing a semi-infinite chain: 



H C = D 



n=0 



where To = 0. The operator f na annihilates an electron in 
a spherical shell centered on the impurity site; this shell 
may be reached, starting from shell 0, by n applications 
of the kinetic energy operator, Eq. ( |2.3a ). The / nCT 's obey 
the anticommutation relations {f} lcn f n '<r} = S n .n'&<y,a' ■ 



The dimensionless coefficients e n and are deter- 
mined by the following recursion relations:E!J 

e„ = (fna\n c /D\f ntr ), (2.15a) 

Tn+ll/n+l.o) = (Hc/D - £ n )\fna) - 7"n|/n-l,cr)> (2.15b) 
1 = Ifn+lAfn+l,*)- (2.15C) 

Here, |/ no ) = /^lO), where |0) is t he va cuum state. 
There is no summation over a in Eqs. ( |2.15| ). 

The first two coefficients generated by the recursion 
relations are 



and 



el— n 

£ = F~ 2 / de e w 2 (e) 



rl—fj, 

n =F- 2 I de (e-£ ) 2 w 2 (e). 



(i+m) 



(2.16) 



(2.17) 



Beyond this point, the expressions for the coefficients 
entering Eq. (2.14) rapidly become complicated. It is 



straightforward to show, however, that if the problem is 
symmetric about the Fermi energy — i.e., if p — and 
w(e) — w{—e) — then e„ = for all n. 

For most functional forms of u>(e), the hopping coeffi- 
cients r n rapidly converge with increasing n to a constant 
value. This prevents faithful approximation of the prob- 
lem using any finite-length chain, because terms in 7i c 
involving sites n which are remote from the impurity are 
just as large as terms involving sites very close to the 
impurity. 



D. Discretization of the conduction band 

Wilson showedEl that, by replacing the continuum of 
conduction band states by a discrete subset, one can in- 
troduce an artificial separation of e nergy scales into the 
hopping coefficients t„ entering Eq. (2.14). This provides 
a convergent approximation to the infinite-chain prob- 
lem using finite- length chains, which correctly reproduces 
the impurity contribution to system properties. Wilson's 
procedure was developed for a flat conduction-band den- 
sity of states, and hence in the notation introduced above, 
for w(e) = wo (a constant). Here, we present a general- 
ization of the method to arbitrary w(e). 

We divide the band into two sets of logarithmic energy 
bins, one each for positive and negative values of e. The 



[£n fLfna + T n (P n Jn-l,, + H.C.)] , (2.14) ^+ 



positive bin (m = 0, 1, 2, . . . ) extends over energies 



l < e < £+, where 



e+ = l-p; e+ = (I-m)A 1 



> 0. (2.18) 



The corresponding negative bin covers the range e m < 
e < e m+1 , where 



-(l+/i); e~=-(l+/i)A 



1—z—m 



m > 0. (2.19) 
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Here A parametrizes the discretization: numerical cal- 
culations are typically performed with A = 2-3, while 
the continuum is recovered in the limit A — > 1. Wilson's 
original treatment of the Kondo problem corresponds to 
setting (i = and = 1. Values z ^ 1 are used in the. 
direct calculation of dynamicaEj and thermodynamic^ 
quantities. (Thermodynamic results are presented in Sec- 
tion VI below.) 

Within the m th positive [negative] bin, we define a 

complete set of destruction operators amir [bmlr] and an 

associated set of orthonormal functions V'am(e) [>pi m (s)], 
q = 0, ±1,±2, all of which vanish for any e outside 
the m th positive [negative] bin. Given such a basis, one 
can write 



oo oo 

<W = £ J2 & { *t(e)alt + ^)b m i] (2.20) 

m— q— — oo 



and 

de e 



tb { - q ' ) *(£)ib {q ' ) (£)a { - q ^a {q ' ) 

Vam V^lVam V fc / "mc u 'ma 



+ ^*(£)4 m \z)&%A- ( 2 - 21 ) 



The key step in generalizing the numerical RG method 
to arbitrary w(e) is the choice of a q — function within 
each bin that has the same energy-dependence as w(s): 



w(e)/F am , 


£ m+l < £ — £ m> 


o, 


otherwise. 


w(e)/F bm , 




o, 


otherwise. 



(2.22a) 
(2.22b) 



The orthonormality condition on these functions implies 
that 



F„ 



dew 2 (e), F{ 



2 

bin 



dew 2 (e). (2.23) 



With this choice, 



.ha = F 



rn— 



,(0) 



'mi",,,,, t Fb m b ma . 



(2.24) 



i.e., the impurity couples only to the q = mode within 
each bin. Following an extension of the reasoning applied 
by Wilson,E3 it can be shown that the coupling be twee n 
modes q ^ q' contained in the kinetic energy [Eq. ( [2.21 )] 
vanishes in the continuum limit A — » 1. To a good ap- 
proximation this coupling can be neglected for A > 1 as 
well. (The "discretization error" arising from this ap- 
proximation is estimated in Section (v|) We therefore 
assume that the q 7^ modes decouple completely from 
the impurity, and contribute to the kinetic energy an un- 
interesting constant term which is dropped henceforth. 
Then, 



H C = DJ2 (Cam a^atl + e bm &<&t 6 W ), (2-25) 
where 



m=0 



-•<„> -F a £ I deem (e), 



^bm Fj 



bm 



de e w 2 (e). 



(2.26a) 
(2.26b) 



Equation fl2.25| ) can now be tridiagonalized using the 
Lanczos recursion relations introduced in the previous 
section. We define 



fner ^ { u nm a ina v nmb m lr) i 



m=0 



where 

UOm — F am /F, VQm — Fbm/F. 

Then Eqs. Q2.15| ) imply that 



' v nm £ bm). 



7"n+l^n+l,m (^6m ^n)^nm ^~n^n — l,m: 
1 = /j( M n+l,m + v n+l,m)- 



(2.27) 



(2.28) 



(2.29a) 

(2.29b) 
(2.29c) 

(2.29d) 



These equations retain the feature of the undiscretized 
conduction band that if fi = and w(e) = w(—e), then 
£„ = for all n. 

As will be discussed in greater detail below, the hop- 
ping coefficients t„ typically decrease like A - ™/ 2 for large 
n, while the on-site energies e n drop off at least this fast. 
For this reason, it is convenient to work with scaled tight- 
binding parameters 



where 



e„ = a- x A n / 2 £ n , t n = a- x k n l\ n , (2.30) 

a= ^(l+A-^Ai- 2 (2.31) 

is a conventional factorilSil which approaches unity in 
the continuum limit A — > 1. With these definitions, the 
discrctized conduction-band Hamiltonian becomes 

oo 

H c = aDY J ^ n/2 [enftjnu 

n=0 



+ t n (flJ n - X , a +R.c.)] 



(2.32) 



In the-soecial case w(e) = wq with fj, = 1 and z = 1, 
Wilsoncj'Ea was able to derive a closed-form algebraic 
expression for the hopping coefficients i„. Bulla et alS3 



6 



have recently presented an ansatz for t n when the scatter- 



ing rate has the pure power-law form given in Eq. (1.3). 
In general, though, the algebraic expressions for the co- 
efhcie nts r apidly become extremely cumbersome, and 
Eqs. ( 2.29| ) must be iterated numerically. A drawback 
of this approach is that .-the recursion relations prove to 
be numerically unstableBtl With double-precision arith- 
metic performed to roughly 16 decimal places, it is typ- 
ically possible to iterate Eqs. ( |2.29| ) only to n = 10 for 
A = 3, and to n = 13 for A = 2. n 

Chen and Jayaprakash have showncl that it is possible 
to reorder the calculation of the t n 's in such a way as 
to circumvent the instability. In this work, however, we 
have adopted a brute-force approach, employing a high- 
precision arithmetic package to compute the coefficients. 
Typically, 120 decimal places suffice for the calculation 
of all coefficients up to n = 30, beyond which point the 
deviation of t n and e n from their asymptotic values is 
insignificant (less than one part in 10 15 ). 



E. The discretized impurity problem 

After discretization of the conduction band, the one- 
impurity Anderson or Kondo Hamiltonian can-be written 
as the limit of a series of finite HamiltoniansBa 



H= lim ak~ N/2 DH 



N — *oo 



N ■ 



(2.33) 



where Hn, describing an (JV+l)-site chain, is defined for 
all N > by the recursion relation 



Hn — A 1//2 -Hat_i + e N f Na f N(X 

+ t N{f Na fN-l,a + H.C.) - Eqn- 



(2.34) 



Here, Eq,n is chosen so that the ground-state energy of 
Hn is zero. The Hamiltonian Ho describes the atomic 
limit of the impurity problem. For the Anderson model, 



Ho = eofl a fo„ + £dn d + Un d] n di 
+ f 1 /2(/t (r d CT+ H.c.)-^G,o, 



where 



aD' 



U = 



u 

aD 1 



r = 







ira 2 D 



while for the Kondo models, 

Ho = (eo + V)flf 0a + jflh** ta >fo* 



s — 



with 



J = F 2 p Q Jo/a, V = F 2 Po Vo/a. 



(2.35) 
(2.36) 

(2.37) 
(2.38) 



In the two-channel variant of the Kondo model, each / 
and p operator in Eqs. (2.34) and (2.37) acquires a chan- 
nel index j, which is summed over. 



An important feature of both the Anderson and 
Kondo Hamiltonians is their behavior under the following 
particle- hole transformations: 



Anderson 
Kondo 



fner * ( 1) fn,(T> 
fna > ( — 1) fn,ai 



-4; 

-st 



(2.39) 



Examination of Eqs. (2.34)^(2.38) indicates that the ef- 
fective values of t n , T, U and J remain unchanged under 
the transformations, but that e„ — ► — e„, e d — > ~(e d +U), 
and V -> -V. 

A symmetric weighting function such that w{—e) = 
w{e) guarantees that e„ = for all n. In this case, one 
sees that the physical properties of the Anderson model 
are identical for impurity energies e d and — (e d + U), all 
other parameters being the same. Thus, it is necessary 
to consider only e d > —U/2 (or e d > —U /2) in ocder to 
fully explore the physical properties of the modelEj 

It should further be noted that, provided w(—e) = 
w(e), the symmetric Anderson model (defined by the 
condition U + 2e d = 0) and the Kondo models with zero 
potential scattering (Vb = 0) are com pletely invariant un- 
der the transformations in Eqs. ( 2.39| ). In cases where the 
impurity scattering rate is regular [w(0) > 0], the pres- 
ence or absence of particle-hole symmetry does little to 
affect the physics. By contrast, this symmetry turns out 
to play a crucial role in determining the strong-coupling 
behavior of systems with power-law scattering. 



F. Iterative solution of the discretized problem 



Th e se quence of Hamiltonians Hn defined by 



Eq. (2.34) can be solved iteratively in the manner de- 
scribed in Refs. ^ and ^3|. The many-body eigenstates 
of iteration N — 1 are used to construct the basis for it- 
eration N. Before each step, the Hamiltonian is rescaled 
by a factor of A 1 / 2 so that the smallest scale in the en- 
ergy spectrum remains of order unity, and at the end of 
the iteration the ground-state energy is subtracted from 
each eigenvalue. This procedure is repeated until the 
eigensolution approaches a fixed point, at which the low- 
lying eigenvalues of Hn are identical to those of Hn+2- 
(The spectra of Hn and Hn+i do not coincide because of 
a fundamental inequivalence between the eigensolutions 
for chains containing odd and even numbers of sites. For 
example, particle-hole symmetry ensures the existence of 
a zero eigenvalue for N+l odd, whereas there is no such 
restriction for N+l even.) 

All the Hamiltonians described in this paper commute 
with the total spin operator 

N 

Stv = fia\cr<j<j> fna* + dl^a^tda', (2.40) 



and with the charge operator 
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Q 



N 



N 

£ 

n=0 



(fLfna ~ 1) + (n d ~ 1). 



(2.41) 



At particle-hole symmetry, Hn commutes with-all three 
components of an "axial charge" operator Jat:E3 



Jz,N ~ ^Qn, 
N 

^ = E(- i r/M;-44> 

n=0 

■7-,iV = (J+,n) ■ 



(2.42a) 
(2.42b) 
(2.42c) 



Thus, the Hamiltonian Hn can be diagonalized inde- 
pendently in subspaces labeled by different values of the 
quantum numbers S z , S, Q, and (at particle- hole symme- 
try) J. Moreover, the energy eigenvalues are independent 
of S z , and when J is a good quantum number they are 
also independent of Q. It is therefore possible to perform 
the numerical RG calculations using a reduced basis con- 
sisting only of states with S z = S and, where appropriate, 
Q = —2 J. (See Refs. ^ and [36] for further details. Note 
that in the two-channel Kondo model, axial charge quan- 
tum numbers and Jw can be defined separately for 
channel j = 1 and 2.) By taking advantage of these sym- 
metries of the Hamiltonian, the numerical effort required 
to diagonalize Hn can be considerably reduced. 

Even with the optimizations described in the previous 
paragraph, after only a few iterations the Hilbert space 
of Hn becomes too large for a complete solution to be 
feasible. Instead, the basis is truncated according to one 
of two possible strategies, which give essentially the same 
results: one either retains the M states of lowest energy, 
where M is a predetermined number; or one retains all 
eigenstates having an energy within some range E c of 
the ground-state energy. A compromise must be made 
between two conflicting goals: accurate reproduction of 
results for the original undiscretized system, favored by 
choosing A to be close to unity to minimize discretization 
error, and by making M or E c large to reduce truncation 
error; and a short computational time, which points to 
large A and small M or E c . Unless otherwise noted, the 
thermodynamic quantities presented in this paper were 
obtained using A = 3 and E c > 25, for which choices 
the primary source of error is the discretizati on. (The 
magnitude of the error is estimated in Section |V C| .) 

In summary, the numerical RG method can be gener- 
alized in a fairly straightforward manner to treat energy- 
dependent scattering of conduction electrons from the 
impurity site. The conduction band is divided into log- 
arithmic energy bins just as for the case of a constant 
scattering rate. However, the mode expansion within 
each energy bin has to be modified in order that the 
impurity couples to a single mode (q = 0). This in turn 
alters the hopping coefficients t n which enter the recur- 
sive definition of the Hamiltonians Hn [see Eq. (2.34)]; 
away from strict particle-hole symmetry, furthermore, 
the hopping terms are complemented by on-site terms of 
the form e n fX a f n <y Finally, a factor F, which depends 



on the overall normalization of the weighting function 
u>(e), is introduced into the starting Hamiltonian Hq [see 
Eqs. (|i3f)-(|Oji)]. 



III. CONDUCTION BAND HAMILTONIANS 

For the remainder of the paper, we focus on systems 
in which the scattering rate T(e) vanishes in power-law 
fashion at the Fermi level. Initially, we consider the sim- 
plest possible form of T(e), given by Eq. (1.3). In the no- 
tation of Section ||, this corresponds to a particle-hole- 
symmetric problem in which fi = and the weighting 
function w(e) satisfies 



w(e) 



|r/2 



0, 



H < i; 

otherwise. 



(3.1) 



Here r can take any non-negative value; r = corre- 
sponds to a constant density of states. Subsequently, we 
shall generalize this form to restrict the power-law vari- 
ation in uu(e) to a finite region, and to allow for particle- 
hole asymmetry. However, the essential phys ics is cap- 



tured by the prototypical function in Eq. (3.1 



In this section, we apply the formalism described in 
Section || to construct and analyze a family of discretized 
conduction-band Hamiltonians based on Eq. ([3.1|), 



N 



= £ A^-")/ 2 [e„/t/ n + t n (/t/ n _ 1 + H.c.)]. 



n—L 



(3.2) 



The parameter L determines the innermost shell ft, onto 
which conduction electrons can hop, so t^ is necessarily 
zero. The case L = represents the free-electron prob- 
lem, while L = 1 and L = 2 will turn out to describe 
electronic excitations at different strong-coupling fixed 
points of the Anderson and Kondo models. In each of 
these limits, spin-up and spin-down electrons decouple 
from one another, so the index a can be dropped. 

All information about the energy-dependence of scat- 
tering from the impurity site enters the discretized An- 
derson and Kondo Hamiltonians through the normaliza- 
tion factor F, and through the tight-binding parameters 
e n and t n . From Eq. (2.5), it is straightforward to see 



that F = 2/(1 + r), while the particle-hole symmetry 



of Eq. (3.1) ensures that e„ = for all n. By contrast, 
the coefficients t n must be determined—numerically, as 
outlined in Section |H D| . (Bulla et alS3 have recently 
deduced an algebraic formula which appears to fit the 
numerical value of t n for all n.) 

Values of t n for r = 0, 0.2, and 1 are illustrated in 
Table Q. Whereas for a constant scattering rate (r = 0), 
t n rapidly approaches unity, in the power-law cases the 
asymptotic value alternates with the parity of n: 



lim t„ 

n — >oo 



t*, 

t*A- r / 2 , 



n odd; 
n even; 



(3.3) 
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scaled hopping coefficient, tn 
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r = 0.2 


r = 1 
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19 


0.9999999991 


1.0174569806 


1.0833333330 


20 


0.9999999997 


0.9115991897 


0.6254627914 


21 


0.9999999999 


1.0174569812 


1.0833333333 


22 


1.0000000000 


0.9115991899 


0.6254627916 


23 


1.0000000000 


1.0174569813 


1.0833333333 


24 


1.0000000000 


0.9115991899 


0.6254627916 


25 


1.0000000000 


1.0174569813 


1.0833333333 



TABLE I. Tight-bi ndin g parameters t n for three different 
powers r entering Eq. (3.1). 



where 



1+r l-A-( 2+r > 



1 + A- 1 2+r l-A-U+f)" 

given by Eq. 



(3.4) 



(3.3) 



The asymptotic behavior of t 
should be contrasted with that obtained for all weighting 
functions in which w(0) is finite and nonzero 



lim e„ = 0, 



lim t n = 1. 



(3.5) 



All such problems_can be showed to exhibit essentially 
the same physics,EEI and in particular to be described 
by the same set of RG fixed points. The functional form 
of w(s) away from the Fermi energy determines only the 
deviations of e„ and t n from their asymptotic values. 
These deviations act as irrelevant perturbations in the 
RG sense. The only exception is eo, which acts as a 
marginal variable, equivalent to an additional potential 
scattering term of the type appearing in Eq. (2.9). The 



scenario of power-law scattering considered in this paper 
is interesting precisely because w(0) — 0, which places 
the impurity problem in a different universality class. 



A. Free-electron Hamiltonian (L — 0) 

The free-electron Hamiltonian, describing the nonin- 
teracting conduction band in the absence of any mag- 
netic impurity level, corresponds to the large-Af limit of 
H { ° } as defined in Eq. |[f). 



Eigenvalues: Numerical diagonalization indicates that 
for large N, the eigenvalues approach limiting values, 
which we denote 



N odd: 
N even: 



j = ±l,±2,.. 
j = 0,±l,±2, 



■ ,±HN+1); 



,±\N. 



Due to the particle-hole symmetry, ?)q = 0. For 
the eigenvalues are well-approximated by 



where 



* = sgn(j) i*A^- 



N odd; 
even. 



(3.6) 



> 1 



(3.7) 



(3.8) 



Eigenvectors: We also consider the single-particle 
eigenstates of , associated with particle operators gj 
(j > 0) and hole operators h_j (j < 0). It will later 
prove useful to have expansions of the original operators 
/„ in terms of these eigenoperators: 



^(JV+l)/2 A 



n.) 



\9i + (-1)"/*}] 



N odd; 



■E£i 2 A y fe + (-l)"/4], Severn 



(3.9) 



It is found, again numerically, that for j ^> 1, the coef- 
ficients Aqj and Aij separate into an A^-dependent pref- 
actor and a part that depends only on the parity of N: 

A nj = A-( 2n+1+r W*a nj , n = Q,l, (3.10) 

where 



[ ao(z )]2 = 1 [! _ A -(l+r) ]A (l+r)(,-l/2) ! 

[ ai (l)] 2 = i[l-A-( 3+r 5]A( 3+r )/ 2 . 



(3.11) 



(The general z-dependence of a\ cannot be written so 
compactly as that of a$.) 

Just as in the standard case r = 0, tho,expansion of all 
other / n 's becomes more complicated£3 For n > 1, f n 
contains components which vary as J^-(. 2m + 1 + r ) N / 4 j for 
m = n, n — 2, n — 4, . . . , nmod2. 



We emphasize that Eqs. (I3/7D, (|3.10|) and (^l|) are 
very good approximations even for comparatively small 
values of Af and j. For A = 3, for instance, these formulae 
hold to at least seven decimal places for all 3 < j < N/2. 
The rate of convergence to the asymptotic forms with 
increasing j seems to be independent of r, at least in the 
range < r < 2. 







B. Symmetric strong-coupling Hamiltonian (L=l) 



Eigenvectors: We can expand the annihilation opera- 



In Section IV E , we shall discuss the symmetric strong- 



coupling fixed point of the Anderson and Kondo models. 
At this fixed point, an infinite coupling between the lo- 
calized level and conduction electrons at the impurity 
site completely suppresses hopping of conduction elec- 
trons onto or off shell 0, i.e., the fo degrees of freedom 
are "frozen out." In this situation, the conduction-band 
excitations of the system are described by an effective 



Hamiltonian H^' obtained by setting L = 1 in Eq. (3.2) 



Eigenvalues: For large A, the eigenvalues of Hjp are 
found to approach limiting values which we denote 



N odd: uj*, j = 0,±1,±2,. 
A even: c2£, j = ±1,±2, . . . 



,±i(iV-l); 



,±±N. 



(3.12) 



Due to the particle- hole symmetry, uj* = 0, while for 
N/2 ^ | j| ^ 1, the eigenvalues are well-approximated by 

u], uj* = sgn(j) (fA- n ' a )AW- v »-\ (3.13) 
where n = min(r, 1). Alternatively, one can write 



u* = a--/ 2 ,;. 



(3.14) 



These expressions point to a significant difference be- 
tween the cases r = and r > 0. In the former instance, 
the strong-coupling energies given by Hjp become iden- 
tical to the weak-coupling energies of ii^-i m ^ ne limit 
of large N. In other words, the transition from weak to 
strong coupling is equivalent to-an interchange between 
the odd- A'' and even- A spectraBa This relation between 
weak and strong coupling no longer holds true in the pres- 
ence of a power-law scattering rate. For r > 1, however, 
an even simpler pattern emerges: uj* — rj* (j ^ 0) and 
= ^-sgn(j) ■ Thus, for superlinear scattering rates, the 
strong-coupling spectrum differs from that at weak cou- 
pling only by the insertion of one additional zero-energy 
eigenstate. 

Phase shifts: The L = 1 eigenvalues can also be ex- 
pressed in terms of s-wave phase shifts So(e) applied to 
the L — spectrum: for A odd, 



uj* = A~ ssn ^' 5o ( Q;A N/2ri ^^r/* , 



(3.15) 



and similarly for A even. The valu es of <5o(e) in th e limits 
e — ► ± can be deduced from Eqs. ( pV7| ) and (|3. 13 ) , while 
the form of the leading corrections away from the Fermi 
energy can be inferred by studying, either pure potential- 
scattering from the impurity sited or power-law mixing 
between conduction electrons and a noninteracting reso- 
nant level (see Appendix |Aj) . The result is 



S (e) = (1 - n) - sgn(-e) 



<D(\e/D\^ 



— r\ \ 



(3.16) 



The main feature of Eq. (3.16) is a jump of (1 — ri)ir in the 
phase shift on crossing the Fermi energy. The interpre- 
tation of this jump will be deferred until Section VD4 



tors f na in terms of single-particle eigenoperators of H N 
Bnogo + E^r 1)/2 B nj [ 9j - (-l) n h% A odd; 



(i). 



fn 



E£i 2 s nj fe-(-i)"/it], 



A even. 
(3.17) 



For sufficiently large N and A/2 ^> j 1, the coeffi- 
cients Bij and Bij separate into an A-dependent pref- 
actor and a part that depends on A only through vn-i- 



f> . - \-(Zn-2+\l-r\)N/A R . 



n = l, 2, 



where 



j3 n - = p n t^ 2n - 2 +\^-r\)(j-VN-i)/2 _ 



(3.18) 



(3.19) 



The parameters /3 n can be determined numerically, but 
we have not obtained algebraic expressions for their de- 
pendenc e on A and r. There is an important exception 
to Eq. ( 3.18 ) for r > 1: in the limit N — > oo, B w for 
A odd \B\\ for A even] approaches a constant value 
which is independent of N . Therefore, the effective scal- 
ing is /i CX A'il-^N/A and h K A (2+|l-r|)iV/4_ Thege 

forms will turn out to have important implications for the 
stability of the symmetric strong-coupling fixed point. 

The approximate expressions for the eigenvalues and 
eigenvectors of do not apply so widely as those for 
L = 0. First, Eqs. ( fUgj ), ( fTlS| ), and ( |3.19| ) are restricted 
not only to j ^> 1, but also to j <C A/2. Second, the rate 
of convergence with increasing A is slower than for L = 
and depends explicitly on r: for instance, the deviation 
of each eigenvalue (uij or ljj) from its large- A limit (uj* 

or uj*) is proportional to A _ l 1_r l Ar / 2 . There is again an 
exception for r > 1: the two smallest even- A eigenvalues 
obey the relation &±i cx A~( r-1 -W 4 , and hence converge 
to their asymptotes (u)±i = 0) even more slowly than the 
other eigenvalues. 

The results above provide clear evidence that a linear 
scattering rate represents a singular case. The expansion 
of /i acquires an A-independent component at r = 1, and 
for this value of r alone several quantities (the eigenvalues 
and eigenvectors as functions of A, the conduction-band 
phase shift as a function of energy) converge in a log- 
arithmic, rather than exponential, fashion. As pointed 
out by Cassanello and Fradkinja this case in some sense 
represents the upper critical dimension of the theory. 



C. Asymmetric strong-coupling Hamiltonian (L — 2) 

It will be shown below that, in most cases, the stable 
strong-coupling fixed point of the Anderson and Kondo 
problems is not described by the conduction-band Hamil- 
tonian Hj^ 1 introduced in the previous subsection. In- 
stead, the system reaches either the frozen-impurity fixed 
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point, described by or the asymmetric strong- 

coupling fixed point, at which both the /q and f\ degrees 
of freedom are frozen out and the low-lying excitations of 

(2) 

the system are described by ' . Due to the asymptotic 
form of the hopping coefficients t n [see Eq. ( |3.3| )1, which 
only depend on whether n is odd or even, the low-energy 
properties of are equivalent to those of H^_ 2 , pro- 
vided that one relabels the operators /„ appropriately, 
i.e., f n at the L = 2 fixed point has the same expansion 
a s fn— 2 at the L — fixed point. The L = 2 eigenstates 
are related to those for L — by a low-energy phase shift 
<5 (e) = Trsgn(-e). 



IV. STRONG- AND WEAK-COUPLING FIXED 
POINTS OF THE IMPURITY MODELS 

In this section, we analyze RG fixed points of the An- 
derson and Kondo m odel s with a power-law scattering 
rate described by Eq. ( 3T ) . The focus is on hosts having 
a single conduction channel, although reference will be 
made in passing to the two-channel Kondo model. 

We consider only those fixed points that can be ob- 
tained by setting each impurity parameter entering Hq 
[see Eqs. (2.35) and ( |2.37 )] either to zero or to infinity. 
In such cases, a local Fermi-liquid description applies: 
the low-energy many-body excitations of the system can 
be constructed as the product of two independent sets of 
single-particle excitations, one set describing the conduc- 
tion band, the other arising from any active impurity de- 
grees of freedom. By studying small deviations from the 
fixed-point Hamiltonian, one can determine the stability 
of the fixed point and the functional dependence of cer- 
tain physical properties. These results can be obtained 
by largely algebraic means. (The only numerical step 



is the derivation of the results presented in Section III 



which involves diagonalization of simple quadratic Hamil- 
tonians.) 



In Section VI, we shall also discuss a number of fixed 



points which appear at intermediate (neither zero nor 
infinite) couplings. Such fixed points are generally non- 
Fermi-liquid in nature, and at present can be studied only 
via a full implementation of the numerical RG scheme 
outlined in Section |l| 



Any deviation from a fixed point Hamiltonian must be 
describable in the form 



SH N = H N - H* = V 70 



7 



7! 



(4.2) 



where 7 is a dimensionless coupling, and 1 is composed 
of operators associated with those degrees of freedom 
(from among f na , d a or s) that remain active at the fixed 
point, multiplied by an overall factor of K N / 2 wh ich re- 
produces the scaling of Hjv implied by Eq. ( 2.34 ). The 
only constraint on the combination of operators entering 
Oj is that the perturbation must preserve all symmetries 
of the original model. 

one can use 



As explained in detail in Refs. 22 and 
the expansion of the operators /„ developed in Section III 
to analyze the stability of the strong- and weak-coupling 
fixed points points with respect to all possible perturba- 
tions. Consider, for example, the weak-coupling limit in 
which the electronic degrees of freedom are described by 



Hff . Equations ( ^9| ) and ( 3. 10| ) imply that the pertur 



bation O v = A N / 2 (f^f a - 5) can be written as A~ rN / 2 
times an A-independent part composed of the single- 
particle and single-hole operators, gj and hj. Making use 
of the effective temperature T/v associated with iteration 
N (see Section [v| for more details) , 



k B T N = aA~ N / 2 D//3, 



(4.3) 



where j3 is a small dimensionless parameter, one sees that 
Oy oc T r , i.e., the perturbation is irrelevant for all r > 0. 

In the remainder of this section, we identify the most 
relevant (or least irrelevant) operators in the vicinity of 
the various Fermi- liquid fixed points. At each fixed point, 
the expansion of f n (n > L) contains a piece which 
varies like 7 1 (™- i )+2 l 1±r l . The dominant perturbations 
are therefore those operators 1 that contain the fewest 
possible fn's, and in which the / n 's that are present have 
the smallest possible indices n. 

We shall present our analysis in the context of the 
nondegenerate Anderson model. Features of the various 
Kondo models will be noted where they are different. 



B. Free-impurity fixed point 



A. Stability of RG fixed points 

Within the nonperturbative RG approach, a fixed- 
point Hamiltonian H* satisfies 



H 



N+2 



H N = H*. 



(4.1) 



In this context, "=" means that two Hamiltonians have 
identical low-energy spectra, and that they share the 
same set of matrix elements of any physically-significant 
operator between their low-lying eigenstates.E3 



The free-impurity or "free-orbital" fixed point of the 
Anderson model corresponds to setting id = U = T = 
in Eq. (2.35). This fixed point, which has no analogue in 
the Kondo models, is described by an effective Hamilto- 



n WC.N — n N 



E, 



G.N ■ 



(4.4) 



Each many-body eigenstate is the product of an eigen- 
state of H $ 
purity level. 



state of and a zero-energy eigenstate of a free hu- 



ll 



By combining the reasoning outlined in the previ- 
ous subse ction with the iV-dependences given in Sec- 
tion [II A, one can identify four operators which are, or 



may be, relevant in the vicinity of the free-impurity fixed 
point: 



Ou 
Or 



k N ' 2 {n d 
k N ' 2 {n d 



1) 
I) 2 



oc T-\ 
oc T~\ 



= A^(/ t < A+H.c. 



T -(l-r)/2 



(4.5) 



= A N / 2 



(fld a + R.c.) oc T-d^)/ 2 . 



Here, £d , Of/, and Or are essentially equivalent to the 
on-site energy, on-site Coulomb repulsion, and hybridiza- 
tion terms (respectively) in the original Hamiltonian, 
while Or d represents correlated hybridization. Of these 
operators, only Ou and Or respect particle-hole symme- 
try and are allowed in the symmetric limit of the Ander- 
son model. Note that O ed and Ou are always relevant, 
whereas Or and Or d are relevant for r < 1 but are ir- 
relevant for r > 1. Since there is at least one relevant 
operator for both the symmetric and asymmetric cases, 
and also for all r, the free-impurity fixed point is always 
unstable. 



C. Valence-fluctuation fixed point 

The valence-fluctuation fixed pointEl of the Anderson 
model corresponds to the original model with e d = T = 0, 
but U = oo. This is clearly not a fixed point of the sym- 
metric Anderson model since U + 2e d ^ 0, and it has 
no analogue in the Kondo models. It is described by 
the same effective Hamiltonian as the free-impurity fixed 
point, but doubly occupied impurity configurations are 
eliminated from the Hilbcrt space. Of the four domi- 
nant perturbation at the free-impurity fixed point [see 
Eqs. ( |4.5| )], only £d and Or survive. Since the former is 
a relevant operator for all r, the valence-fluctuation fixed 
point is always unstable. 



None of the perturbations in Eqs. (4.5) is allowed at 
the local-moment fixed point. Instead, the dominant per- 
turbations are as follows: 



Oj = K N/2 fl \cr aa , f 0a , ■ s cx T r , 

Oy = A N / 2 (flf 0a 1) ex T, 

O tl = K N l 2 (flf la +U.c.) ex T 1+r , 

Uo = A»/ 2 (flf „-l) 2 ex T 1+2r 



(4.7) 



Here, Oj and Oy describe exchange (Kondo) scattering 
and pure potential scattering, respectively; Ot t is a term 
from the kinetic energy; and Ou represents a Coulomb 
interaction between two conduction electrons in shell 0. 
Only Oy breaks particlc-holc symmetry. 

In the standard case r = 0, exchange and pnj.e-. 
potential scattering are marginal; further analysisE3'E3 
reveals that antifcrromagnetic [ferromagnetic] exchange 
is marginally relevant [marginally irrelevant], hence the 
fixed point is unstable [stable]. For r > 0, by contrast, 
all perturbations are irrelevant, so the local-moment fixed 
point is stable irrespective of the sign of J. This is the 
first of several important differences between the fixed- 
point behaviors for r = and r > 0. 

We note that for the two-chan nel Kondo model, each 
of the operators 7 listed in Eq. (4.7) should be replaced 



by a pair of operators, Oy = O^ ± 7 2 ', where Oy 
(j = 1, 2) is identical to 7 except that all its / oper- 
ators carry a channel label j. Throughout this paper it 
is assumed that the two conduction channels couple to 
the impurity spin with equal strength, in which case only 
the symmetric operator 7 + ^ can enter SHjy. In addition, 
one can construct allowed perturbations that contain /'s 
belonging to both channels. At the local-moment fixed 
point, the leading perturbation of this type is 

O v , = A^ 2 (/t lCT / 0lCT - 1)(/oWo2. - 1) oc T 1+2r . (4.8) 

Similar remarks concerning the two-channel Kondo 
model apply at each of the remaining fixed points de- 
scribed in this Section. 



D. Local- moment fixed point 

The local-moment fixed point corresponds to the orig- 
inal Anderson model with T = and U = —2e d = oo. 
The eff ectiv e Hamiltonian at the fixed point is still given 
by Eq. (4.4), but only singly occupied impurity states are 
allowed, so there is a decoupled spin-one-half degree of 
freedom, 



4 



2 " jUjU 7 "O/z' ? 



(4.6) 



localized at the impurity site. This is the weak-coupling 
fixed point of the Kondo mode ls, i.e., it corresponds to 
setting J = V = in Eq. fl2.37|). 



E. Symmetric strong-coupling fixed point 

The symmetric strong- coupling fixed point is obtained 
by setting T = oo in Eq. fl2.35| ) while keeping s d and 
U finite, or by setting J = oo in Eq. ( 2.37 ). We first 
consider the ground state of the atomic Hamiltonian Hq. 
In the nondegenerate Anderson model and the conven- 
tional s = i Kondo model, any moment at the impurity 
site is competely screened by /o electrons. (In the An- 
derson model, the ground-state impurity occupancy (n d ) 
varies continuously with the parameters e d and U. For 
any given set of couplings, however, there is a unique 
ground state, and hence no residual impurity degree of 
freedom.) The s = 1 Kondo model and the two-channel 
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Kondo model have spin-one- half ground states^ in the 
former instance, the impurity is underscreened, in the 
latter, it is overscreened. 

In all cases, an infinite gap separates the ground 
state(s) from all other eigenstates of the atomic Hamil- 
tonian. As a consequence, the foe degrees of freedom are 
frozen out and the effective Hamiltonian becomes 



H 



(i) 

N 



Eg. 



jv- 



(4.9) 



Based on the results of Section III B| , the dominant 
perturbations at this fixed point are 



-ri 



-ri 



0. h = A N / 2 fl \cToo- ha' ' T OC T 

Vl = h N ' 2 {flh„ 1) oc T- 
O u% = k N '\flf la -lf oc T 1 
O t2 = A w / 2 (/ 1 t J 2ff +H.c.) oc r (i-r 1+ |i-r|)/ 2i 

(4.10) 

where, as before, r\ = min(r, 1). Oj t , a Kondo- like oper- 
ator involving the residual spin t, is present only in the 
underscreened and overscreened models. Oy 1 describes 
nonlocal potential scattering of electrons in shell 1 from 
the impurity site. Both Oj 1 and Ov 1 are relevant pertur- 
bations for all r > 0. Ojj 1 , representing Coulomb repul- 
sion between f\ electrons, is relevant for r > i. Finally, 
the kinetic energy term Ot 2 is marginal for r = 1, but is 
irrelevant otherwise. 

The fixed point is stable for r > if, and only if, 
three conditions are satisfied: the impurity moment is 
exactly screened (to rule out Oj 1 as an allowed pertur- 
bation); the power r is less than i (to ensure that Ojj 1 is 
irrelevant); and the problem exhibits particle-hole sym- 
metry (so that Ovi is disallowed). Thus, one sees that 
the symmetric strong- coupling fixed point is generically 
unstable. This represents another significant departure 
from the standard case r = 0, in which the fixed point 
is always marginally stable, except in.oyerscreened prob- 
lems, where it is marginally unstablcBd 



F. Asymmetric strong-coupling fixed point 

The asymmetric strong- coupling fixed point of the An- 
derson model and of the single-channel Kondo models is 
obtained in the same fashion as the symmetric strong- 
coupling fixed point considered above, with a further 
cond ition : either the coefficient e± entering Eqs. ( 2.32| ) 
and (2.34) is made infinite; or the model Hamiltonian Hn 
is augmented by a term ViOvn where Oy 1 is defined in 
Eqs. (4.10) and |Vi| — ► oo. As a result, the f\ degrees of 
freedom are frozen, in addition to those associated with 
shell and the impurity. The effective Hamiltonian be- 
comes 



The two fixed points described by ei = +oo and 
e\ = — oo (or by V\ = +oo and V\ = — o o) ha ve different 



ground-state charges as defined in Eq. (2.41), but they 
are otherwise physically equivalent and will henceforth 
be treated as a single fixed point. 
The dominant perturbations are 



Oj 2 

o U2 



= A^/ 2 /U<wiW 

= A^HfLha+H-c) 
= A N/2 (fLh.-l) 2 



oc T r , 
oc T r , 

rpX-\~r 



(4.12) 



x 



oc 



where f describes a residual spin-one-half degree of free- 
dom present only in the underscreened Kondo model. 
Since these operators are all irrelevant, the asymmetric 
strong- coupling fixed point is stable for all r > 0. 

In the standard case (r = 0) , the symmetric and asym- 
metric strong-coupling fixed points represent two points 
on a continuous line of marginally stable fixed points de- 
scribed by a family of effective Hamiltonians 



H* N (Vi) = H$ + Vxk^tflh,, -l)-E t 



G.N- 



(4.13) 



These fixed points share essentially the same physical 
properties, independent of the value of V\. The effect 
of a power-law scattering rate to destroy all but two of 
these fixed points and to make the case V\ = unstable 
with respect to the breaking of particle-hole symmetry. 

The two-channel Kondo model also has a stable, 
strong-coupling fixed point at J = oo, e\ = ±oo. How- 
ever, we have not found any choice of the bare parame- 
ters Jo and Vq that produces flow to this limit, in which 
the ground state carries a residual spin-one-half degree 
of freedom. Instead, it is helpful to consider the Hamil- 
tonian Hn 



see Eqs. (4.1C) and the 



last paragraph of Section IV D]. The asymmetric strong- 
coupling fixed point of interest is reached by first setting 
J — oo and V = to lock the impurity into an over- 
screened spin doublet, and then taking the simultane- 
ous limits J\ — > oo and V\ — > ±oo in such a way that 
1/2 < \V\\/J\ < 3/4. Under this prescription, the impu- 
rity combines with shells and 1 to produce two degener- 
ate ground-state configurations carrying quantum num- 
bers (S,QW,Q( 2 )) = (0,Tl,0) and (0,0, Tl)- This pair 
of spinless states represents a "flavor-one-half" degree of 
freedom. [The generators of electron flavor symmetry 
are obtained from the standard spin generators by inter- 
changing spin and channel indices: "f<-> 1, J,<-> 2. For ex- 
ample, the z-componcnt of flavor measures the difference 
between the number of electrons in channels 1 and 2, i.e., 

jfiavor = Uq(1) _ q(2))_] fixed pQint j g described by 



the effective Hamiltonian H\ sc N defined in Eq. (4.11) 



the leading irrelevant perturbations are and a flavor 
analogue of . 



H*asc,n — Et N 



(2) 



E, 



G,N- 



(4.11) 
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G. Frozen-impurity fixed point 

The frozen-impurity fixed point of the Ande rson model 
is obtained by setting id = +00 in Eq. ( [2.35 ). Here, the 
impurity level becomes completely depopulated, an d th e 
excitations of the system are just described by Eq. (4.4). 
The leading perturbati ons at this fixed point are Oy, 
O tl and Ou from Eqs. (4/7), so the fixed point is clearly 
stable for all r > 0. 

The electronic excitations at the frozen-impurity and 
asymmetric strong-coupling fixed points are described 
by a nd respectively. As pointed out in Sec- 

tion III C , these two Hamiltonians can be made equiva- 
lent by a suitable relabeling of the operators /„. The 
leading irrelevant perturbations about the two fixed 
points become identical under this relabeling. Thus, 
the frozen- impurity and asymmetric strong- coupling fixed 
points are physically equivalent, up to a shift in the 
ground-state charge. In treating the Anderson Hamilto- 
nian, it will prove more convenient to refer to the frozen- 
impurity fixed point (since (nd) can be made arbitrar- 
ily small in this model), whereas the asymmetric strong- 
coupling fixed point more naturally describes the Kondo 
models (which correspond to the limit rid — 2s). 



V. THERMODYNAMIC PROPERTIES 

This section is primarily concerned with the numerical 
and analytical calculation of the contribution made by 
magnetic impurities to various thermodynamic proper- 
ties. First, though, we remark briefly on the thermody- 
namics of the pure host Fermi systems. 



A. Host Thermodynamic Properties 

As we have emphasized in the Introduction, the im- 
purity properties of the models we consider depend on 
the conduction-band density of states and the energy- 
dependent hybridization only in the particular combina- 
tion r(e) = irp(e)t 2 (e). However, in order to compute the 
properties of the pure system in the absence of magnetic 
impurities, it is necessary to specify p(e) explicitly. If 
the power-law energy-dependence of the scattering rate 
arises solely from the hybridization, then the host proper- 
ties will be those of a conventional metal. Here we focus 
on the opposite limit, more appropriate for describing the 
gapless systems listed in the Introduction, in which the 
hybridization is essentially const ant and the density of 
states has the form given in Eq. (El). In this case, unit 
normalization of p(e) implies that 

1 + r 



Pa = 



2D 



(5.1) 



With these assumptions, it is straightforward to show 
that for k B T <C D, the host entropy, specific heat capac- 
ity, and static susceptibility are given by 



SM (knT 



£7(0) 



k B 

k B T X {0) No 
2 



2N (l + r) (2+r)0(l+r) 



l+r 



k B T 
D 



l+r 



(l + r) 



knT 



l+r 



(5.2a) 
(5.2b) 
(5.2c) 



(gp B ) 2 2 ^" 1 ■' V D 
Nq is the number of unit cells making up the solid, and 
4>(x)=( 1 (x + l), $(x)=x[C 1 (x)-C2{x)], (5.3) 
where, for all x > and all positive integers n, we dcfincEl 
f°° u 1 " 1 

Ux) = L du WTrr- (5 - 4) 

(The functions <j) and <f> will also enter the impurity prop- 



erties calculated in Section VD 



One sees from Eqs. ( |5.2| ) that the exponent r which 
determines the density of states is directly reflected in the 
temperature-dependence of the host pxoperties. For later 
reference, we define the host Wilsono (or Sommerfeld) 
ratio, 



R 



(o) 
w 



4tt 2 fc|T X (0) 



(l + r) 



t5o 3 (gp B ) 2 C^ 3(l+r)(2+r)0(l + r)' 

(5.5) 



Since (f)(1) = 7r 2 /12 and (f)(1) = |, Eq. ( |575|) reduces in 
the limit r — > to the standard result, R' 



(0) 

w 



B. Impurity Thermodynamic Properties 

The impurity contribution to a thermodynamic prop- 
erty A is defined to be the change in the total measured 
value of A brought about by adding a single impurity 
to the system. Each such contribution can be computed 
from an expression of the form 



A imp = (A) imp = (A) - (A) 

'Tr(Ae- /l » ait ) -Tr (Ae' 13 " 11 ^ )] , (5.6) 



= lim 



where A is an operator which depends on the property 
of interest, 



13 N = a DK- N l 2 /(k B T) 



(5.7) 



is the natural energy scale of iteration divided by the 
thermal energy scale, and "Tr " means a trace taken over 
an impurity-free system. 

For example, the impurity contributions to the entropy 
and the specific heat are obtained as 



q V-L lm p 

Jimp — QJ^' 



d 2 F 

a- - t 1 



imp 



(5.8) 
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Here, Fj mp is the difference between the total Helmholtz 
free energy of the system with and without the impurity: 

-Fimn = -knT In Zw 



with 



= lim fcsTln ( Z^P/Zn 



(5.9) 



(5.10) 



Another quantity of interest is the impurity contribu- 
tion to the zero-field magnetic susceptibility, given by 



imp / n O I ry\ 

= \O z l £j /imp 



lim 



Tr(S 2 



(5.11) 



Tr n f Sle-f 3 ""* 



7 (o) 

N 

where /is is the Bohr magneton, g is the Lande g factor 
(assumed to be the same for conduction and localized 
electrons), and S z is the z-component of the total spin of 
the system. The quantity SkgTximp equals the square of 
the effective moment contributed by the impurity to the 
system. 

The numerical RG formulation provides a controlled 
approximation for computing the impurity contribu tion s 
to each thermodynamic property according to Eq. ( |5.6| ) . 
The method does not yield reliable results for (A) or 
{A)o separately, even though these are the values that 
would be have to be measured experimentally in order to 
determine Aj mp . 



Numerical Evaluation of Impurity 
Thermodynamic Properties 



In Section VI 



we present thermodynamic properties 
obtai ned v ia the direct numerical evaluation of Eqs. ( |5.9| ) 
and ( |5.11 ). This subsection briefly reviews some of the 
technical details of these calculations. 

The general strategy for computing thermodynamic 
properties using the discretized Hamiltonians Hn is as 
follows: One first selects a value for the dimensionless 
parameter (3. Then for each iteration N = 0, 1, . . . , one 



assigns the re sult of Eq. ( 5.E ) to the temperature defined 
through Eq. ( |5.7|) by the condition (3n = (3- This gives 
the quantity (Ah mp at a sequence of temperatures Tjv 
satisfying Eq. (|4.3|). The T/v's are equally spaced at inter- 
vals of h In A on a logarithmic scale. If desired, this "grid" 
of temperatures can be refined by using several different 
values of j3 at each iteration. The choice (3 = (3~oA^^ 2M 
for j = 0, . . . M (we have used M = 4) proves convenient 
because the corresponding temperature grid {Tnj} con- 
tains the redundancies Xjv,o = T^+i.m- The discrepancy 



between the two independent evaluations of a thermo- 
dynamic property at the same temperature provides a 
useful measure of the error in the result. 

Since the small est energy scale of Hn is of order unity, 
one expects Eq. (5.6) to provide increasingly reliable re- 
sults for A > 1 as (3 becomes much smaller than unity. 
However, there is another factor which militates against 
taking f3 <C 1. Limitations of computer time and memory 
permit the retention only of those states having an energy 
within E c of the ground state. In order to minimize the 
contribution of the missing states to {A)i mp , one wants 
(3E C to be as large as possible. In practice, therefore, 
(3 is chosen as a compromise to take a value somewhat 
smaller than one. The results presented in Section VI 
were calculated for A = 3 or 9, retaining all eigenstates 
up to a dimensionless energy 25 and using values of (3 
between 0.6 and 0.6A -1 / 2 . It is shown in Ref. [2^ now 
one can calculate corrections to compensate for such rel- 
atively large values of /3. Our studies indicate that while 
the corrections are formally of order (3/ A ~ 0.1-0.2, they 
have small prefactors which reduce the overall shift in 
Ci mp and Ximp to less than one part in 10 3 . Since this 
level of error is smaller than that arising from the dis- 
cretization of the conduction band, we have neglected 
the /3/A corrections. 

As a practical matter, Eqs. (|5.8|) are not used directly 
to evaluate Simp and Ci mp . A more accurate evaluation 
of the entropy, which avoids numerical differentiation, 
exploits the relation 



Sh 



k B (((3H/Zy iinp - 0F\mp) ■ 



(5.12) 



It is likewise possible to obtain the specific heat without 
differentiation, through the equation 



Cm 



(j3H/Zy) -[({3H/Z) imp 

imp 



(5.13) 



but the results turn out to be rather prone to discretiza- 
tion error. All plots of the specific heat presented below 
were instead obtained using C; mp = TdS- lulv / dT with a 
simple two-point approximation to the derivative. 

As mentioned above, the thermodynamic quantities 
presented in this paper were obtained using A = 3 or 
A = 9 and an energy cutoff E c > 25. With these choices, 
the primary source of error is the discretization. One 
of the main effects of working with a value of A greater 
than unity is the introduction into (^4)i mp of oscillations 
which are periodic in InT. The oscillations have a period 
In A and a magnitude proportional to exp(— 7r 2 / In A). 
Oliveira and Oliveira have shownEj that these oscillations 
can be greatly reduced by averaging values of (A)i mp 
computed for different band discretization parameters z 
(see Section [II Dj ). We have employed four z'a (0.5, 0.75, 
1, 1.25) in obtaining the results presented below. 

Another consequence of the band discretization is a 
reduction in the effective coupling between impurity and 
delocalized degrees of freedoms Study of a discretized, 
resonant-level model with a power-law scattering ratec2l 
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indicates that the most faithful descr ipti on of the c on- 
tinuum problems described by Eqs. (2/7) and ( 2.11 ) is 
obtained by premultiplying the parameters To, po^O) and 
PqVq entering the discretized calculations by a factor 



A(A,r) 



■l_A-(2+r)" 


l+r 




l + r 


2 + r 




1 


-A-( 1+r ) 



2+r 



In A, (5.14) 



l + C(lnA) 2 



for A 1. 



The remainder of this section summarizes properties 
of the five distinct fixed points discussed in Section IV . 
Perturbative corrections to the fixed-point properties are 
presented for each of the three stable (or conditionally 
stable) regimes. Certain technical details have been rel- 
egated to Appendix IbI 



1. Free-impurity fixed point 



For r = 0, Eq. ( |5.14 ) reduces to the standard result given 
in Ref. EM In the remainder of this paper, we quote the 
continuum equivalent of each coupling. Thus, any nu- 
merical data labeled with a particular value of r , poJo, 
or poVo were actually comp uted by substi tutin g ATq, 
Ap J , or Ap V Q into Eqs. fl2.36| ) or Eqs. fl2.38| ). Note 
that parameters desc ribin g the impurity alone, i.e., 
and U entering Eqs. ( [2.36 ), do not have to be corrected. 

The measures outlined in the preceding paragraphs 
greatly reduce, but cannot completely eliminate, dis- 
cretization errors in the computed thermodynamic prop- 
erties. We estimate on the basis of limited calculations 
performed for other values of A that for A = 3, the over- 
all error in the impurity contributions to the suscepti- 
bility Tx or the entropy S is less than 5%, while that 
for the specific heat C is less than 10%. It should be 
emphasized that these are errors in absolute quantities 
at finite temperatures. Zero-temperature properties and 
exponents describing ratios of properties at different tem- 
peratures or couplings generally have much smaller errors 
(below 1%). Indeed, fixed-point properties can be com- 
puted to better than 1% using values of A as large as 
10, with a considerable reduction in the numerical effort 
compared to that required for A = 3. 



D. Perturbative Evaluation of Impurity 
Thermodynamic Properties 

In the vicinity of any of the fixed points described in 
Section IV , perturbation theory can be applied to the ap- 



propriate effective Hamiltonian to obtain analytical ex- 
pressions for thermodynamic quantities as functions of 
the couplings 7 which parametrize the deviation from 
the fixed point [see Eq. (4.2)]. Once perturbative expres- 



sions have been obtained for the discretized version of 
the problem (A > 1), they can be extrapolated to the 
continuum limit (A = 1). 

A similar perturbative treatment of the standard 
Kondo and Anderson models is described in detail in 
Refs. ^2] and ^3], respectively. The extension to sys- 
tems with a power-law scattering rate is conceptually 
straightforward but algebraically laborious. One novel 
feature is that in certain physical regimes the dominant 
temperature-dependences derive from second-order cor- 
rections to the fixed-point properties, whereas in the 
standard case (r = 0) is is not necessary to go beyond 
first order in perturbation theory. 



The free-impurity fixed point describes a decoupled or- 
bital which has four degenerate states. The impurity 
contributes to the thermodynamic properties as follows: 



Simp = k B In 4, 



Cimp 0; 



^bTx 



imp 



(.9mb y 



|. (5-15) 



2. Valence-fluctuation fixed point 

The valence-fluctuation fixed point describes a decou- 
pled orbital which has three degenerate states, double 
occupation being forbidden. The impurity properties in- 
clude 



S n 



fee m 3, Ci m p = 0, 



^bTx 



imp 



(.9MB ) 2 



i (5.16) 



3. Local-moment regime 



One can apply standard perturbative methodsE3£3 to 
the effective Hamiltonian — H1^ c jv+S 7 7^7 > where 
the fixed-point Hamiltonian H* wc N is given by Eq. ( |4.4| ) 
and the perturbations 1 are those defined in Eq. ( [4.7] ). 
For the nondegenerate Anderson model or the exactly 
screened Kondo model one finds, to lowest order in each 
of the perturbative couplings, that 



* = In 2 - 8hl3 N A-^ N / 2 £ a 0j a ljPj 



k B T 

j 

+ 2(y 2 + Aj 2 )feA--(^M (5.17) 
1 \''> / i/fc ''>■ J 



and 



imp 
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(5.18) 



.Simp = k B hi 2, Cj 



-n k BT Ximp _ 1 fco'Vl 
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The corresponding results for the underscreened model 
are 



S'imp — In 3, C; m p — 0, 



&4 (5 - M) 



a 0jPjPj 



Pj-Pj 



P N (l-6pjPj) 



In all four models, the leading corrections at low tem- 
peratures take the form 



Here. 



A*Si m p, ACi m p 



X 



T ri+r , A(T Ximp ) oc T r (5.25) 



Pj 



1 



Pi 



1 



(5.19) 



is the occupation probability of a fermionic state having 
energy r]*, and the indices j and k run over the range 1 
to (N+l)/2, inclusive. We have omitted the lowest order 
contribution of Ojj to -Fi mp because it is highly irrelevant 
(oc T 3+4r ). The equations above are written for N odd; 
similar expressions hold for N even. 

In the continuum lim it (N — ► oo , A — > 1, and /3jv <C 1), 
the sums entering Eqs. ( 5.17|)_a nd (5.18) can be evaluated 
algebraically (see Appendix |]|) . The resulting equations, 
valid for all positive r except r = 1, are 



for all positive r ^ 1. 

For the special case r = 1, the second-order terms in 
-Pimp and Ximp acquire logarithmic corrections, necessi- 
tating the replacement 

ip(r) 



In 



D 



Pk B T 



(5.26) 



± imp 



= In 2 - 4tiV(l + r)(3 + r)0(l+r) 



k B T 



D 



l+r 



+ 2(y 2 + ^ J 2 ){l+r) 2 ^(r) 0(n+r) 



ri+r 



(5.20) 



in Eqs. (5.2C) and ( 5.211). Here, (i is the small param- 
eter introduced in Eq. ( |4.3| ), the precise value of which 
cannot be determined uniquely within the present for- 
malism. (See Appendix [b| for further discussion.) As a 
result, there is some uncertainty in the thermodynamic 
properties, but it is clear that ASimp and ACi mp must 
contain contributions proportional to T 2 log T and oth- 
ers varying like T 2 . The ratio of the prefactors of these 
terms will determine whether or not the logarithmic cor- 
rection is observable at temperatures of physi cal interest. 
We shall return to this point in Section VI A . 



and 
(9^b) 2 
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V 2 

+ yfl+fll'W 0(ri+r) 



k B T 
D 



ri+r 



(5.21) 



Both <f>{x) and <f>(x), defined in Eqs. (5.3), vary smoothly 
with x and are of order unity over the range < x < 3. 
By contrast, the function 



4- Symmetric strong-coupling regime 

The methods of the previous section can be applied to 
compute impurity properties at the symmetric strong- 
coupling fixed point, plus the leading corrections for 
those cases in which the fixed point is stable, i.e., for 
the particle-hole-symmetric Anderson model and the ex- 
actly screened Kondo model, both with r < \. Here, 
first-order perturbation theory suffices, yielding the ex- 
pressions 



In 4 + 4 In q t - 4 \npj 

l 3 

- M 2 p N k-^ N ' 2 Pilfcm (5.27) 



k B T 



7T rix 
, — tan — 

ijj(r) = { 2 2 
(r-1)- 1 



for < r < 1, 
for r > 1, 



(5.22) 



and 

k B T\i_ 



has a simple pole at r = 1. 

Similar calculations can be performed for the s = 1 
and two-channel Kondo models. To summarize the re- 
sults, the fixed-point impurity properties of the Ander- 
son model and of the two s — | Kondo models are just 
those expected for a decoupled spin-i impurity: 



(9P>b ) 



E _ E 



PjPj 



2t 2 p N K-^ N ' 2 Y,PuP2iqiqi{qi ~ ffl) 
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where pj and pj are given by Eq. (5.19), and qi and qi 
are the analogous quantities defined for the eigenenergies 



°f Hssc,n : 



91 



1 



9t = 



1 + e -^ Nu l 



(5.29) 



The index I runs from 1 to (N— 1)/2, inclusive, while j 
still runs from 1 to (iV + l)/2. 



The sums entering Eqs. (5.27) and (5.28) can be per- 
formed algebraically for values of A close to unity. Ex- 
trapolation of the resulting expressions (see Appendix |b|) 
to the continuum limit gives 
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— n In 4 — ^ib\bi x 
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and 
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D 
D 



1— 2r 



where 



6 2 (r) = Jim 



2^ 



A->1 (2n + l-r) In A 
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(5.32) 



It is found numerically that b\ and 62 are of order unity, 
at least for all < r < i, the range over which the 
symmetric strong-coupling fixed point is stable. 

Thus, the fixed-point impurity properties of the An- 
derson and exactly screened Kondo models are 



S imp = 2r 1 k B In 2, 



k B Tx 



imp 



(.9Ms) 2 8 ' 
those of the underscreened Kondo model are 



S imp = (l+2n)k B In 2, 



k B Txi: 



2+ri 
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and those of the overscreened model are 

k B Txim P _ 1+n 



(l+4n) k B In 2, 
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(5.34) 



(5.35) 



In all four models, C\ mp = at the fixed point, and the 
leading corrections to the fixed-point properties vary as 



AS'impj ACi mp cx T 



l-r 



A(Tximp) oc T 



l-2r 



(5.36) 



The fixed-point properties and temperature exponents 
above agree with-,those obtained for < r < 1 by Chen 
and JayaprakashEl fox. the exactly screened Kondo model 
and by Bulla at aZ.Ej for the Anderson model. [At ex- 
tremely low temperatures (T/D < IO -10 ), the latter au- 
thors identify a T 1_2r variation in the specific heat. Such 



a ter m can arise only from the operator Ojj 1 defined in 
Eqs. ( |4.10| ). Under conditions of strict particle-hole sym- 
metry- however, 0\j 1 cannot contribute to the specific 
heat.EJ We suspect, therefore, that the T 1_2r behavior 
is a numerical artefact.] In addition, Bulla et al. find 
that for < r < 1 the impurity spectral function A(uj) 



varies like 



in the limit 



0, in contrast to the 



Lorentzian form found in the standard case r = 0. 

At the symmetric strong-coupling fixed point, one 
sees that each conduction band makes a contribution to 
the entropy and to the susceptibility equal to r\ times 
that of an isolated level described by Eq. (2.2b) with 
Ed = U = 0. Furthermore, the impurity dens ity o f states, 
Pimp(e) = Tr~ 1 d5o/de computed using Eq. ( 3.16 ), has a 
delta-function peak of weight r\ at e = (see also Ref . |j]) . 
These observations suggest the phenomenological inter- 
pretation that a fraction 7*1 of a conduction electron from 
each band occupies a decoupled level of zero energy, the 
remaining fraction 1 — r± presumably being absorbed into 
a many-body resonance centered on the Fermi energy. 

The properties described above appear to be highly 
anomalous. It should be emphasized, though, that since 
the impurity has no internal degree of freedom at the 
fixed point (at least in the Anderson and exactly screened 
Kondo models), it acts only to exclude conduction elec- 
trons from its immediate vicinity. The fixed point behav- 
iors are simply those of independent electrons subjected 
to a phase shift. Indeed, the noninteracting (U = 0) limit 
of the Anders on m odel reproduces the low-energy phase 
shifts of Eq. ( 3.16 ) and hence yields pre cisely the ther- 
modynamic properties described in Eqs. ( |5.33| ). Further- 
more, the noninteracting model is shown in Appendix ^ 
to exhibit a spectral function A(oj) oc |jj|l 1 ~' r l~ 1 1 in agree- 
ment with the U > results of Ref. [L0|. We conclude 
that the symmetric strong- coupling fixed point embodies 
a natural generalization of standard Fermi-liquid physics 
to gapless hosts. 



5. Asymmetric strong-coupling/frozen-impurity regime 

The thermodynamic properties at the frozen-impurity 
fixed point of the Anderson model and at the asymmetric 
strong-coupling fixed point of the exactly screened Kondo 
model are simply 



Cimp 0, 



(5.37) 



The corresponding properties of the underscreened 
Kond o mo del are those of a free spin-one-half, given by 
Eqs. (5.23). The overscreened Kondo model has a decou- 
pled flavor-one-half degree of freedom, which contributes 



Su 



In 2, C h 



0, T X , 



0. 



(5.38) 



The corrections to the fixed-point values can be ob- 
tained from the corrections at the local-moment fixed 
point by the replacements J — > 0, V — > V2, ti — > £3, and 
Uq — * U2- The results are of the form 
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AS imp , ACimp oc T ri+r , A(Tximp) oc T r > +r . (5.39) 

Note that these expressions do not extrapolate to r = 0, 
in which limit the leading corrections are linear in T. 

For r = 1, the thermodynamic properties above should 
be supplemented by terms proportional to T 2 logT. In 
this same case, Cassanello and FradkinQ have found loga- 
rithmic corrections to the Kondo temperature, the static 
susceptibility, and the T-matrix. These authors point 
out that under such circumstances, multiple energy scales 
enter the problem and physical quantities are no longer 
controlled by the Kondo scale alone. 

The results of the previous two paragraphs imply that 
for the Anderson and exactly screened Kondo models, 
both Ci m p and Tximp approach zero as T T1+r . This is 
the only instance among all the fixed points and models 
discussed in this paper in which a system described by a 
scattering exponent r > exhibits a nontrivial impurity 
Wilson ratio, 



R 



w 



^B^Ximp 



lim — 



(5.40) 



Examination of Eqs. ( |5.20 ) and ( |5.21 ) shows that for 
r > 1, Rw is a function of V$ 1^3 as well as of r, and 
is thus expected to depend on the bare couplings (Fo, 
ed, and U for the Anderson model; J and V for the 
Kondo model). For < r < 1, however, the leading 
contribution to both the susceptibility and the specific 
heat is proportional to V 2 2 alone, so the Wilson ratio takes 
a universal value, 



Rw(r) 



n 2 (/>{2r) 



6r(l + 2r)(t){2r)' 



(5.41) 



This expression will be compared with the host Wilson 
ratio (r) in the next section. 

VI. NUMERICAL RESULTS 

This section presents numerical RG results obtained 
using the formalism described in the earlier parts of this 
paper. We concentrate primarily on pure power-law scat- 
tering rates of the form of Eq. (1.3). At the end of the 
section, we discuss the effect of various modifications to 
the scattering rate, including the introduction of particle- 
hole asymmetry and the restriction of the power-law vari- 
ation to a pseudogap region around the Fermi energy. 

For simplicity, we shall henceforth set ks = gfJ-B = 1- 
(We remind the reader that the g-factor is assumed to be 
the same for localized and conduction electrons.) 

A. Anderson model 

This subs ecti on treats the nondegenerate Anderson 
model, Eq. (2.1), restricted to the domain U > 0. For 



pure power-law scattering rates of the form of Eq. ( |l.3| ) 
it is nec essary to consider only > —{7/2 (see Sec- 
tion HE). Most of the results will be presented for the 
extreme cases U = —2td and U = oo which, respectively, 
preserve and maximally break particle-hole symmetry. 

We first examine in some detail the properties of the 
model for a fixed exponent, r — 0.2. In particular, we 
show how the variation of the thermodynamic properties 
with decreasing temperature can be interpreted in terms 
of crossovers between various of the fixed-point regimes 
enumerated in Section IV. This interpretation can in 



some instances be corroborated by computing the rela- 
tive populations of the four impurity configurations. We 
then explore some of the systematic changes that take 
place as r is varied, focusing on the progressive damp- 
ing of impurity charge fluctuations and the consequent 
suppression of the Kondo effect. Finally, we relate the 
preceding results to the simple scaling theory of Ref. 0. 



1. Pure power-law scattering rate with r — 0.2 

Figures |l] and || provide an idea of the range of possible 
behaviors in the temperature-variation of the impurity 
magnetic susceptibility and specific heat. In these plots, 
U and Tq are held fixed, and each curve represents a dif- 
ferent value of ed- (To prevent overcrowding, we gener- 
ally place a symbol at only one in every six temperature 
points along each curve when plotting thermodynamic 
quantities. The line connecting the symbols results from 
a fit through the complete data set.) 

The case U — oo is shown in Fig. [|. At very high 
temperatures, T D (not shown), the properties are 
close to th ose of the valence-fluctuation fixed point (Sec- 
IVC|): Txi,„p ~ 1/6 and C ilnp w 0. Once T drops 



tion 

below D, the properties become sensitive to the posi- 
tion of the impurity level relative to the Fermi level. If 
Ed is positive or weakly negative (e.g., see the curve for 
(d/ D = —0.1), Tximp falls monotonically with decreasing 
temperature, indicating a crossover from valence fluctu- 
ation to the frozen- impurity regime (Section IV G ). This 
crossover is accompanied by a peak in the specific heat 
representing a loss of impurity entropy equal to In 3. 

For more negative values of ed, Tximp initially rises 
as T falls below D, but at lower temperatures it drops 
back towards zero (see the curves for €d/D = —0.15 
and -0.175 in Fig. |). The rise can be associated with 
a crossover from valence fluctuation to local-moment be- 



havior (Section IV D| ) , even though T^imp does not climb 



all the way to 1/4, the value characterizing a free spin 
s = i. The subsequent drop in Txi mp signals a sec- 
ond crossover to the frozen-impurity fixed point as the 
impurity becomes Kondo-screened. The specific heat 
shows two well-defined peaks, corresponding to the two- 
stage quenching of the impurity entropy from In 3 at 
the valence- fluctuation fixed point to In 2 in the local- 
moment regime to zero at strong coupling. This double- 
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FIG. 1. Impurity susceptibility Tximp and specific heat 
Cimp vs temperature T for the infinite-C/ Anderson model 
with values of as labeled and a pure power-law scattering 
rate specified by r = 0.2 and To = 0.1 D. T he legend shows 
the local-moment fraction /lm [see Eq. (pTlt)] for each curve. 



peak structure may be taken as a signature of the Kondo 
effect, just as it is in a system with a flat scattering rate. 

If the impurity level lies far below the Fermi energy 
{td/D = —0.3 and —0.5 in Fig. |l|), then as the temper- 
ature falls, Tximp rises monotonically towards the local- 
moment value of 1/4. The crossover between from va- 
lence fluctuation to local-moment behavior is marked by 
a single peak in the specific heat as the impurity loses 
the entropy associated with the empty-impurity configu- 
ration. In contrast to the standard case r = 0, where the 
system eventually flows to strong coupling for any choice 
of bare impurity parameters, these curves show that for 
r = 0.2 and ed sufficiently negative, an unscreened spin 
survives on the impurity site down to absolute zero. This 
is made possible by the local stability of the local-moment 
fixed point for all r > (see Sectio nJIV D ) . 

The final curve shown in Fig. ^fJeJ/D = -0.223) 
quickly reaches a plateau at Tximp ~ 0.2, and remains 
there down to at least T = 10 -6 D. This behavior is not 
comp atib le with any of the fixed points described in Sec- 
tions IV and [y|. Moreover, it is achieved only by careful 
tuning of td for given r, U and Tq. We therefore inter- 



pret it as evidence for an unstable, intermediate-coupling 
fixed point — the manifestation in the Anderson model 
of the fixed point identified by Withoff and Fradkincl in 
the Kondo model. This identification is discussed further 
in Section 
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As mentioned in the introduction to this subsection, 
it can also be useful to examine the ground-state impu- 
rity configuration. At temperatures sufficiently high that 
the system is in the free-impurity (T > U) or valence- 
fluctuation (T < U) regime, one expects the Anderson 
model to exhibit charge fluctuations between subspaces 
labeled by different values of n d - However, once T drops 
below the effective values of both r and |e<j|, real charge 
transfer is frozen out and (at least for ed > —U/2), 
the local level should be well-described by one of three 
configurations: (i) a local moment (only states with 
rid = 1 are significantly populated); (ii) an empty im- 
purity ((nd) ~ 0); (iii) a mixed-valence configuration (in- 
volving significant occupation of states with more than 
one nd value). 

The three low-temperature impurity configurations 
can be differentiated using the local-moment fraction, 



/lm = Km(n d 



2n d -\n d i}. 



(6.1) 



For U = oo, /lm is identical to the ground-state im- 
purity occupancy (nd). For U = —2ed, however, (nd) al- 
ways equals one (due to particle- hole symmetry), whereas 
/lm varies from zero to one-half to one as is changed 
from +oo to to — oo. 

The legend on Fig. [I] lists the local-moment frac- 
tion for each of the curves. The value /lm = 0.63 
for ed/D = —0.1 places the impurity within the mixed- 
valence range. All the remaining curves have /lm > 0.75, 
signaling the existence of a well-developed local mo- 
ment. It is interesting to compare this information with 
that provided by the thermodynamic properties. Note 
that the cases e d /D = —0.10 and —0.175 result in flow 
to the same fixed point, despite having the different 
ground-state impurity configurations. Conversely, the 
local-moment configuration present for e d /D — —0.175, 
—0.223, and —0.30 is nonetheless associated with three 
distinct fixed-point behaviors. These observations serve 
to emphasize that /lm, which probes only the local im- 
purity configuration in the ground state of the system, 
is complementary to the fixed-point analysis, which is 
based on the low-energy excitations above the many- 
body ground state. 

The susceptibility plots for a symmetric Anderson im- 
purity in Fig. [2] exhibit many of the features found for 
U = 00. In particular, large, negative values of 
drive the system to the local-moment regime, and there 
exists a critical value (ed/D = —0.40) for which T^imp 
remains approximately equal to 0.2 down to very low 
temperatures. However, the susceptibility curves for 
fd/D > —0.40 head towards a low-temperature limit of 
TXimp ~ 0.025 instead of zero. This is precisely the be- 
havior expected at the symmetric strong-coupling fixed 
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10" 4 10" 3 
k B T/D 

FIG. 2. Impurity susceptibility Tximp and specific heat 
Cimp vs temperature T for the symmetric Anderson model 
with values of as labeled and a pure power-law scattering 
rate specified by r = 0.2 and To = 0.1D. 



point (see Sections [VE| and |VD4 ) , which exhibits a sus- 
ceptibility Tximp = ri/8, where n = min(r, 1). This 
departure from the case U = oo demonstrates the im- 
portance of particle-hole (a)symmetry in determining the 
low-temperature properties. 

The specific heat curves in Fig. ^ are qualitatively sim- 
ilar to those in Fig. [j], a double peak again indicating the 
formation and subsequent Kondo-screening of a local mo- 
ment. However, differences in quantitative features such 
as the area under each peak reflect the fact that the cases 
U — oo and U — —2e c i have different fixed points in both 
the high- and low-temperature limits. 

The approach to the low-temperature fixed points 
is highlighted in Fig. [|, which shows the deviations 
of Tximp and Ci mp from their zero-temperature lim- 
its. Examples are provided of both the infinite-?/ and 
symmetric cases, and for flow both to strong coupling 
(ed/D = —0.05) and to the local-moment fixed point 
{ed/D = —0.50, —0.65). (A symbol is placed at every 
second data point in this figure.) 

Each curve in Fig. [| is labeled with an exponent ob- 
tained by fitting the low-temperature data to a power law 
in T . The exponents for the U = oo curves are all close to 




k B T/D 

FIG. 3. Deviation of properties from their fixed-point 
values: impurity susceptibility A(Txim P ) and specific heat 
ACi m p vs temperature T for the Anderson model with a pure 
power-law scattering rate specified by r = 0.2 and To = 0.1-D. 
Straight lines represent fits to the data points over three 
decades of temperature. Each curve is labeled with its fit- 
ted slope (estimated error ±0.01). 



r or 2r, in good quantitative agreement with Eqs. ( 5.25 ) 
and (5.39). For a symmetric impurity there are g reate r 
depa rture s from the asymptotic forms in Eqs. ( 5.25j ) 
and j5 36D , The exponent of 0.64 ± 0.01 for A(Tximp) 
in the symmetric strong-coupling regime (e^ = — 0.10Z3) 
probably reflects the admixture of a substantial resid- 
ual T 1 ~ r contribution into the leading T 1 ~ 2r term [see 
Eq. ( 5.31| )]. The deviation from T 2r behavior in ACi mp 
near the local-moment fixed point (e<j = — 0.65-D) can 
also be attributed to incomplete convergence; fits limited 
to the lowest decade of temperatures for which reliable 
data are available yield an exponent of 0.39 ± 0.02, com- 
pletely consistent with perturbation theory. 

Figure [| shows that the low-temperature state can also 
be tuned by varying U at fixed Tq and e^. Particularly in- 
teresting are the three curves that show signs of entry to 
the local-moment regime and subsequent Kondo screen- 
ing. The middle curve (U/D = 0.5) displays the anoma- 
lous properties associated with the symmetric strong- 
coupling fixed point (Tximp = fi/8, S- lmp — r\ In 4). The 
curves on either side (U/D = 0.4 and 0.6) eventually en- 
ter the frozen-impurity regime (Tximp = 0, S- lmp = 0), 
but only after lingering close to the symmetric fixed point 
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FIG. 4. Impurity susceptibility Tximp, specific heat Ci mp , 
and entropy Si mp vs temperature T for the Anderson model 
with td = — 0.25.D, values of U as labeled, and a pure 
power-law scattering rate specified by r — 0.2 and To = 0.1 IX 



over some range of temperatures. (The effect is especially 
pronounced for U/D = 0.6. This curve has only just be- 
gun to move away from the symmetric fixed point at the 
lowest temperatures shown.) 



2. Trends with increasing r 

In the paragraphs that follow, we examine some of the 
systematic trends that arise when r is varied. We begin 
with the position of the unstable intermediate-coupling 
fixed point which separates the stable weak- and strong- 
coupling basins of attraction. Discussion of the physical 
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FIG. 5. Position of the intermediate-coupling fixed point, 
F c vs r for the Anderson model with a pure power-law scat- 
tering rate. Solid lines are provided as a guide to the eye. 



properties of this fixed point will be deferred until Sec- 
tion 
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As shown in Figs. the intermediate fixed point can 
be located by adjusting e ( i at fixed Tq and U, or by tun- 
ing U while holding Tq and e<j constant. The third pos- 
sibility is to define a critical scattering rate r c , such that 
for Tq > r c the system flows to strong coupling, while 
for Tq < T c the low-temperature physics is governed by 
the local-moment fixed point. Figure || plots F c /|ed| as 
a function of r for two fixed values of ed, and for both 
U = oo and U = —2ed- The two infinite-?/ curves coin- 
cide almost perfectly and are roughly linear in r. For a 
symmetric impurity, the dependence on r and e<j is more 
complicated, the most notable feature being the diver- 
gence of the critical scattering rate as r — > |, beyond 
which point the strong-coupling fixed point is completely 
inaccessible. These trends will be discussed further in 



Section VI B 3 



We now turn to the strong-coupling behavior governed 
by the frozen- impurity fixed point. The perturbation the- 
ory of Section |v| indicates that in this regime the impurity 
Wilson ratio Rw defined in Eq. (5.40) takes a universal 
(although r-dependent) value over the range < r < 1. 
The Wilson ratio can be obtained numerically from the 
computed values of Txi mp and Ci mp . To the accuracy 
that we can achieve, our results for 0.2 < r < 0.8 confirm 
that there is indeed a single value of Rw for each r. 

Figure ^ compares the best value of Rw (i.e., the value 
with the smallest estimated error) determined numeri- 
cally using a discretization parameter A = 3 with the 
continuum perturbative value, Eq. (5.41). The two sets 
of results agree to within 1.5% for r < 0.5, 5% fo r r = 0.6 
and 8% for r = 0.8. (As shown in Section VD 5, the sub- 



leading corrections to T^imp and Ci mp are smaller than 
the leading terms by a factor proportional to T 1_r . Thus, 
as r increases Rw must be calculated at progressively 
lower temperatures, producing larger errors.) Fig. ^ also 
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FIG. 6. Wilson ratio Rw vs r for the infinite-C/ Anderson 
model with a power-law scattering rate. Individual points , 
representing the impurity Wilson ratio defined in Eq. (5.40), 
were determined numerically using a discretization parame- 
ter A = 3. The solid line shows the perturbative result of 
Eq. ( |5,4l| ) for the continuum limit, A — > 1. The dashed line 
plots the host Wilson ratio in the absence of impurities, cal- 
culated assuming a power-law density of states [see Eq. (fTq)]. 
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FIG. 7. Boundaries of the mixed-valence regime, e* d+ 
(open symbols, dashed lines) and e* d _ (filled symbols, solid 
lines) vs To for the mfinite-(7 Anderson model with several dif- 
ferent pure power-law scattering rates. Individual points were 
determined by the criteria /lm = 0.75, 0.25 while the line s 



were obtained using poor-man's scaling (see Section VI A 3) 
Dashed lines from top to bottom correspond to r = 0, 0.25 
0.5 and 0.75. For r > 1, the scaling theory predicts e* d+ = e a _ 



plots the host Wilson ratio R w \ defined in Eqs. (5.5) 



If the deviation from unity of the ratio Rw I Rw can ^ e 
taken as a measure of the impurity-induced many-body 
effects, it appears that these effects weaken or even van- 
ish as r approaches 1. 

As noted above, the strong-coupling regime can be 
divided — based on the ground-state impurity configu- 
ration — into empty-impurity, mixed-valence and local- 
moment subregions. Figure 3 of Ref. ^ shows the vari- 
ation of the impurity occupancy with td for U = oo, 
To = 0.1D and for values of r ranging from to 2. In 
all cases, (n^) (= /lm for U = oo) increases from zero 
to one as the impurity level moves from far below to far 
above the Fermi energy. The effect of a power-law scat- 
tering rate is to narrow the range of td over which the 
crossover takes place from an empty impurity, through 
mixed valence, to a full local moment. 

In order to quantify the narrowing of the mixed- valence 
regime, we define e* d+ and e* d _ as the values of the im- 
purity energy which result in a ground-state occupation 
/lm = 0.25 and 0.75, respectively. These energies, which 
we take to represent upper and lower bounds on the 
mixed- valence regime, are plotted in Fig. [?]. It is clear 
that the mixed- valence region of parameter space shrinks 
monotonically as the power r increases at fixed Tq. 

The local-moment fraction for a symmetric impurity 
with Tq = 0.1D is plotted as a function of the impurity 
energy in Fig. |^. (Only negative values of td are shown. 
Due to particle- hole symmetry, /lm is mapped to 1 — /lm 
when the sign of td is reversed.) Just as for the infinite-C/ 
case, the range of Cd over which /lm takes values between 
0.75 and 0.25 decreases dramatically with increasing r. 



Figure H plots /lm — \ against — on a log-log scale. 
For r < 0.5, /lm clearly dips downward to approach the 
value \ linearly as — * 0. The curves for r = 0.75 
and 0.9 also show signs of the same behavior, although 
the linear regime is pushed to much smaller \e<i\- By 
contrast, the curvature of the data for r = 1 and r = 2 
suggests that /lm approaches a value greater than ^ as 
Ed —>■ 0~, and hence undergoes a discontinuous jump 
when the impurity level passes through the Fermi energy. 

Figures show that as the exponent r describing the 
power-law scattering rate increases, there is a progres- 
sive shrinking of the mixed-valence region of parameter 
space in favor of the local-moment and empty-impurity 
regimes. This trend is a natural consequence of the de- 
pression of the scattering rate near the Fermi level, which 
strongly inhibits mixing between the impurity level and 
low-energy conduction states. 

The numerics also reveal a second trend with increas- 
ing r, namely the progressive disappearance of the Kondo 
effect. For r > 0.5 it proves almost impossible to find 
any region of parameter space within which the ther- 
modynamic properties indicate that entry to the local- 
moment regime is followed by a crossover to either of the 
two strong-coupling fixed points. This is illustrated in 
Fig. |l0|, which shows the impurity susceptibility and spe- 
cific heat for a linear scattering rate. As noted above, 
the case r = 1 is of particular interest because it may de= 
scribe a magnetic impurity in a ci-wave superconductord'cl 
and in flux phases of two-dimensional electrons. 

Consider first the solid curves in Fig. [l(] representing 
the case U = oo. For ed/D > —0.037, the impurity level 
is almost unoccupied in the ground state (/lm ^ 1/4), 
and the system crosses directly from valence fluctuation 
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FIG. 8. Local-moment fraction /lm vs impurity energy 
for the symmetric Anderson model with pure power-law scat- 
tering rates specified by To = 0.1 D and several different values 
of r. The lines are provided as a guide to the eye. 
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FIG. 9. Same as Fig. |], plotted on a log-log scale. 



to the frozen-impurity regime; there is neither a peak in 
TXimp nor a double peak in d mp to signal the Kondo ef- 
fect. For ed/D < —0.040, the renormalization is from va- 
lence fluctuation to the stable local-moment fixed point, 
and the T = ground state has an unquenched spin 
at the impurity site. Between these behaviors lies an 
unstable intermediate-coupling fixed point, obtained by 
tuning the impurity level to ed/D w —0.0379. The ab- 
solutely flat Tximp curve in this case indicates that the 
fixed point is reached directly from the high-temperature 
regime, rather than from the local-moment regime as was 
found for r — 0.2 (see Figs. |] and ||). 

Figure [w] also provides one representative example of 
the properties exhibited by a symmetric impurity with 
a linear scattering rate (dashed line). As far as we have 
been able to determine, any negative value of ed, however 
small in magnitude, results in flow to the local-moment 
fixed-point. This observation is consistent both with the 
absence of any finite critical coupling r c for r > A (as 
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FIG. 10. Impurity susceptibility Tximp and specific heat 
Ci m p vs temperature T for the Anderson model with a linear 
scattering rate specified by r — 1 and To = 0.1D. Data points 
connected by solid lines correspond to U = oo and values of 
Ed shown in the legend. The dashed line shows the symmetric 
case U/2 = -e d = 0.0AD. 



shown in Fig. and with the evidence that there is a 
jump in the local-moment fraction at ed — for r > 1 
(see the discussion of Fig. |9|) . 

Figure [ll] illustrates the approach of the r = 1 thermo- 
dynamic properties to their zero-temperature values. As 
discussed in Section V D , a linear scattering rate admits 
logarithmic corrections to power-law behaviors at both 
the local-moment and frozen-impurity fixed points. For 
T < 10 -3 .D, ATximp is linear in T on appro ach t o the 
local-moment fixed point [as predicted by Eq. (5.21)] and 
quadratic (with no significant component of T z log T) 
in the frozen-impurity regime. The specific heat is ex- 
pected to behave as T 2 logT at both fixed points. The 
deviations from a pure T 2 form are only weak for the 
frozen-impurity example (ed = —0.0115) but are consid- 
erably stronger for the two flows to the local-moment 
fixed point. In all three cases, the lowest decade of tem- 
perature before the data become too noisy for a reliable 
fit seems to be consistent with a T 2 log T form, although 
higher-accuracy calculations would be needed to rule out 
completely any other behavior. 
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k B T/D 

FIG. 11. Deviation of properties from their fixed-point 
values for the Anderson model with a linear scattering rate 
specified by r — 1 and To = 0.1D. Top: Impurity suscepti- 
bility A(Tximp) vs temperature T. Straight lines represent 
fits to the data over the range k B T/D < 1(T 3 . Each curve is 
labeled with its fitted slope (estimated error ±0.01). Bottom: 
Specific heat ACi mp /T 2 vs T, with straight-line fits to the 
data over selected temperature ranges. 



3. Comparison with scaling theory 




FIG. 12. Scaling of the effective impurity energy td 
with the temperature T, shown schematically for (a) r = 0; 
(b) < r < 1. Renormalization of Id from its bare 
value td begins on entry into the valence-fluctuation (VF) 
regime, at D ~ Ud = min(!7, D). Scaling ends at a 
crossover to local-moment (LM), empty-impurity (EI) or 
mixed-valence (MV) behavior. A power-law scattering rate 
flattens the trajectories €d(T) and also moves the boundary 
of the MV regime to the left, reducing the range of bare im- 
purity energies (e^_ < ed < e d+ ) which yield MV behavior. 



The numerical results presented above can be com- 
pared with the predictions of an approximate analytical 
treatment of the problem based on the poor-man's scak 
ing technique first developed for the Kondo problem.EZI 
Jeffcrsonc3 and HaldaneCJ applied this method to the 
nondegenerate Anderson model with a flat scattering rate 
(r = 0). The main effect of many-body interactions was 
found to be a temperature-dependent shift in the effec- 
tive energy of a nonsymmetric impurity level from its 
bare value ed to 

e d (T) = e d +~m(^), r = 0. (6.2) 

Here Ud = min(J7, D) is the energy scale below which 
many-body effects come into play. Neither U nor T is 
significantly renormalized. (Note that particle-hole sym- 
metry prevents renormalization of ed for a symmetric im- 
purity. In this case, the upper bound Ud on the range 



of T within which renormalization occurs essentially co- 
incides with the lower bound |e<j|.) 

In this scaling picture, real charge fluctuations on the 
impurity site are expected to become frozen out around 
a temperature Tp = max(|e<j|, T). (As usual, we assume 
that U > 2\e d \.) If T F = -e d > T, then only the singly 
occupied configurations will be significantly populated, 
and the impurity will possess a local moment, while an 
empty impurity will result if Tp — +Ed 3> T. Finally, if 
Tp = T >• | e^l, then the ground state will have mixed 
valence. Figure |l2|(a) provides a schematic representa- 
tion of the renormalization of ed in the case r = and 
illustrates the crossover from valence fluctuation into the 
three low-temperature regimes. 

The poor-man's scaling treatment- was recently ex- 
tended to a power-law scattering rate.cl The depression of 
the scattering rate at low energies can be represented as 
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a renormalization of the parameter Tq entering Eq. (1.3): 

f (T) = T • (T/D) r , T < D. (6.3) 

This in turn feeds back to produce a smaller renormal- 
ization of e d than occurs in the case r = 0: 



e d (T) = e d + ^ 
nr 



Ud 
D 



r > 0. (6.4) 



The crossover temperature characterizing the freezing- 
out of real charge fluctuations on the impurity site must 
be redefined to be 



T F =max(|e d (T F )|,r (T F )). 



(6.5) 



The scaling of e d (T) for r > is shown schematically 
in Fig. |l^(b). Note that the trajectories are flatter than 
in (a), reflecting the reduction in the many-body shift 
of e d . Moreover, the energy-dependence of the scattering 
rate pushes the vertical line T — Tq(T) significantly to 
the left, thereby shrinking the range of e d values which 
result in a mixed-valence ground state. (For r > 1 the 
vertical line is driven to T = 0, and the mixed-valence 
regime disappears altogether.) 

Bounds on the mixed-valence region. Within the poor- 
man's scaling approach, upper and lower bounds on the 
mixed- valence regime, e d+ and e* d _, respectively, can be 
defined as the bare values of e d for which the scaling tra- 
jectories pass through the upper and lower intersections 
between dashed lines in Fig. [l2|, i.e., as roots of the im- 
plicit equation 



T (\e d (T F )\)=±a ± e d (T F ), 



(6.6) 



where a + and a_ ar e po sitive constants of order unity. 
The solutions of Eq. (6J3) are 



D V 7rr 



r 

a ± D 



l/(l-r) 



To 
■nrD 



Ud 
D 



(6.7) 



Figure superimposes these bounds for the mixed- 
valence region in the case U = o o (plotted as lines) 
on those defined in Section VI A 1 based on the com- 



puted value of the local-moment fraction /lm (individ- 
ual symbols). The choice a_ = 0.9 brings the alterna- 
tive definitions of e* d _ into good quantitative agreement. 
The dashed lines representing e* d+ in Fig. |?], computed 
for a+ = 0.6, fit the numerical data reasonably well for 
small values of To but deviate considerably for stronger 
impurity scattering. This discrepancy is not especially 
surprising, given the approximations inherent to poor- 
man's scaling and the degree of arbitrariness present in 
both definitions of e* d+ . 

For the symmetric case shown in Figs. || and|^, the cri- 
terion /lm = 0.75 for the border of the local-moment 
regime can be compared with the scaling definition 
f (|e2-l) = — a - e d-- (Recall that the level energy does 
not renormalize at particle- hole symmetry, so e d = e d .) 



A value of a_ « 0.4 seems to yield reasonable agreement 
with the numerics, at least for r < 0.5. 

Exchange on entry to the local-moment regime. The 
scaling theory also sheds light on the disappearance of 
the Kondo effect with increasing r. On entry to the local- 
moment regime, the Anderson model can be mapped 
onto the Kondo problem by projecti ng int o the subspacg 
in which n d = 1. Applying Eqs. (|2.10 ), one obtainsO 
a dimensionless exchange having a power-law energy- 
dependence, 



^(eM?)J(e,e')=p J \ee'/D 2 \ r/2 , 



where 



PoJo 



2f (T F ) fj_ 



U -T f 



(6.8) 



(6.9) 



In this instance, Eq. (|6.5| ) reduces to T F = —e d (T F ). 

For given impurity parameters (e d , U , and To), the ex- 
change decreases with increasing r, as shown in Fig. 2 
of Ref. ||. Both the depression of T(T F ) and the weaker 
renormalization of e d (which increases T F ) contribute to 
this effect. A conservative bound on the multiplicative re- 
duction factor for poJoi obtained by neglecting the renor- 
malization of e d altogether, is \e d /D\ r . This sharp reduc- 
tion of pqJq with increasing r militates strongly against 
any Kondo |effect, because (as first shown by Withoff 
and Fradkincl), an impurity spin becomes screened only 
if Jo > J c , where po J c ~ r. 

In summary, the poor-man's scaling analysis captures 
the essential features of the Anderson problem with a 
power-law scattering rate, and provides a convenient the- 
oretical framework for understanding the numerical re- 
sults. A number of quantitative features, however, are 
not accounted for correctly within the scaling approach. 



B. Screened Kondo model 

The conventional s = i ("exactly screened") Kondo 
model with a power-law exchange coupling [Eq. ( |6.g| )l 
has been discussed extensively in Refs. and reviewed 
briefly rin the Introduction above. (Some of the papers 
citedHa formally treat the degenerate Anderson model 
using large-A^ methods; however, the focus throughout 
is the Kondo physics of the local-moment regime.) We 
recall that the novel feature of the model is the existence 
of a finite exchange coupling J c {r, Vq) which separates a 
region of parameter space within which the impurity spin 
becomes asymptotically free (Jo < J c ) from another in 
which the impurity moment is quenched. The latter case 
is governed by two distinct fixed points: under conditions 
of strict particle-hole symmetry,tj the low-temperature 
susceptibility is a universal function of T/Tk, where 
Tk oc I Jo — Jcl 1 ^; otherwise, the properties are deter- 
mined by two independent energy scales u 

Rather than attempting to provide a comprehensive 
treatment of the exactly screened Kondo model, which 
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FIG. 13. Low-lying eigenvalues E* at the parti- 

cle-hole-symmetric intermediate-coupling fixed point of the 
exactly screened Kondo model, plotted vs the exponent r de- 
scribing the power-law exchange. Data computed with A = 3 
are shown both for N even (open symbols) and for N odd 
(filled symbols). Solid lines are provided as a guide to the 
eye. The curves are extrapolated to the weak-coupling values 
at r — and to symmetric strong coupling at r = 0.5. 



would necessarily duplicate much previously published 
work, we focus in this section on three topics that have 
received little attention. First, we address the proper- 
ties of the Kondo model at the critical coupling. We 
show that over a range of r there in fact exist two dis- 
tinct intermediate-coupling fixed points — one accessible 
only under conditions of exact particle-hole symmetry, 
the other reached when this symmetry is broken. Second, 
we present results for a local response function which is 
a candidate order parameter for the critical behavior at 
the intermediate-coupling fixed point (s). Finally, we ex- 
amine the relationship between the Kondo and Anderson 
models in systems with a power-law scattering rate. Al- 
though each RG fixed point of the Kondo model is equiv- 
alent to a fixed point of the Anderson model, we argue 
that two models are independent to a greater extent than 
is the case for r = 0. 



1. Properties of the intermediate- coupling fixed point 

Unlike the weak- and strong-coupling fixed points dis 
cussed in Section [V 



the intermediate-coupling fixed 
points of the Anderson and Kondo Hamiltonians are not 
amenable to conventional perturbative methods, and our 
knowledge of their properties comes entirely from numer- 
ical studies. This subsection presents numerical RG re- 
sults for the Kondo model obtained by tuning the ex- 
change parameter Jo at fixed r and Vq to lie very close to 
the border between the strong- and weak-coupling basins 
of attraction. (These results will be compa red with those 
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FIG. 14. Impurity susceptibility Tximp and entropy 
Simp at the intermediate-coupling fixed points of the exactly 
screened Kondo model, plotted vs r. Solid lines are pro- 
vided as a guide to the eye. The curves are extrapolated 
to weak-coupling at r = 0, and for Vo — are extrap olate d 
to the symmetric strong-coupling values [given by Eq. (5.33), 



dashed lines] at r = |. Inset: Approach of Tximp to its value 
at r = i, plotted on a log- log scale. 



Since we focus on zero-temperature (fixed-point) prop- 
erties, it is not necessary to average over different dis- 
cretizations of the conduction band (as described in Sec- 
tion VC). The results discussed below were obtained 
using a sin gle di scret ization, corresponding to z = 1 in 
" and (|2~19| ). 



Eqs. (gd 



for the Anderson model in Section |VIB 3[ ) 



We begin by considering the particle-hole-symmetric 
problem corresponding to Vb = 0, for which the critical 
coupling J c (r) is-figite only over the range < r < r max , 
where r max = i .ETlI When the bare exchange Jo is tuned 
precisely to J c (r), the problem flows to the intermediate- 
coupling fixed point mentioned previously. Figure [l3| 
plots the lowest-lying eigenenergies at this fixed point 
over the range 0.05 < r < 0.48, while Fig. |l| shows 
(filled symbols) the impurity contribution to the total 
susceptibility and entropy for 0.1 < r < 0.49. As r is 
varied from zero to \ each curve interpolates smoothly 
between the corresponding weak-coupling and strong- 
coupling values. Over the entire range of r, the com- 
puted entropy at the intermediate fixed point remains 
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within 1% of In 2. The inset to Fig. 14 shows that the 
approach of Tximp to its value at r = ¥is described by a 
power law: Tximp - j$ oc (± - r) v with v = 0.54 ± 0.05. 
All the quantities plotted in Figs. [l3] and [l4| vary linearly 
at small r, and are consistent with the divergence of the 
critical coupling J c as r — ► r max = |. 

Away from particle-hole symmetry, a finite critical cou- 
pling J'(r, Vq) can be identified for any r > 0. For r < i, 
J' deviates smoothly from its particle-hole-symmetric 
value J c (r) as Vo is increased from zero. The initial slope 
dJ'/dVo can be of either sign, but for strong potential 
scattering (which disfavors the presence of a single con- 
duction electron at the impurity site) , J' invariably rises 
sharply (see Fig. 1 of Ref. 0). 

The nature of the fixed point reached for Jo = J' 
is found to be fundamentally different for small and 
large values of r. Particlc-holc symmetry proves to be 
marginally irrelevant for 0<r<r*«0.4, over which 
range the intermediate fixed point is identical to that ob- 
tained for Vq = 0. For r > r* , by contrast, all positive 
bare values of Vb result in flow to a new, particlc-hole- 
asymmetric fixed point located at Vo = V c , Jo = J' c = 
J'(r,V c ). Thus, over the limited range r* < r < r max , 
there exist two distinct intermediate- coupling fixed points, 
located at J c and J' c . Potential scattering is marginally 
relevant at the former, whereas the coupling Vo — V c is 
marginally irrelevant at the latter. (Yet another fixed 
point is reached for all Vo < 0, but since it is trivially 
related to that for Vq > by particle-hole exchange, we 
treat this pair as being physically equivalent.) 

Support for the statements contained in the previous 
paragraph comes from the impurity susceptibility and en- 
tropy shown in Fig. [l4|. Properties computed for a fixed 
value of the potential scattering, p Vo = 0.1, are essen- 
tially indistinguishable from their Vo = counterparts 
for all r < 0.4. Beyond this point, T^imp and S'imp at the 
particle-hole-asymmetric fixed point rise monotonically 
with increasing r. The properties for poVb = 0.5 (not 
shown in Fig. |l4|) are identical to those for poVo = 0.1, 
at least to within our estimated accuracy. 

The many-body spectrum provides additional evidence 
for the existence of dual fixed points over the range 
r* < r < r max . Certain pairs of charge-conjugate states 
which are necessarily degenerate for Vo = are split at 
the asymmetric fixed point. Figure [T| plots the split- 
ting AE* of the lowest pair of affected states, both for 
odd- and even-numbered iterations N, over the range 
0.38 < r < 1. The vanishing of AE* defines the critical 
value r*. For r < 0.45, we find that AE* oc (r - r*) u \ 
with r* = 0.375 ± 0.002 and v* = 0.67 ± 0.15. It seems 
probable that the magnitude of the splitting is directly 
proportional to the critical potential scattering, at least 
for small V c , but we have no proof of this conjecture. 

Figures 16(a), (b) and (c) summarize the effect 
of particle-hole asymmetry on the exactly screened 
Kondo model for a fixed value of r, where < r < r* , 
r* < r < r max , and r > r max , respectively. These figures 
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FIG. 15. Energy splittings AE* of the lowest pair of 
charge-conjugate eigenstates at the particle-hole-asymmetric 
intermediate-coupling fixed point of the exactly screened 
Kondo model, plotted vs the reduced exponent r/r* — 1, for 
A = 3 and r* = 0.3754. Straight lines show fits to the left- 
most four data points for odd- and even-numbered iterations 
of the numerical RG method. 



sketch the RG flow of the effective couplings Jo and Vb 
on the plane Vi = for fixed r. [Here V\ measures the 
strength of potential scattering exper ienced by electrons 
in the Wilson shell /i; see Eq. ( |4.10| ). Note that asym- 
metric strong coupling corresponds to |Vi| = oo. Out-of- 
plane flows towards this fixed point are represented by 
broken lines.] As r increases from zero, the intermediate- 
coupling fixed point moves steadily to the right along the 
horizontal axis. At r — r* , the particle-hole-asymmetric 
fixed point separates from that for Vo = 0. The two fixed 
points grow further apart as r rises towards r max = -|, 
at which point the symmetric fixed point merges into 
the strong-coupling limit. Beyond r max , the asymmet- 
ric fixed point remains at finite couplings. We believe 
that it continues to move upward and to the right with 
increasing r. 



2. Local impurity susceptibility 

In a recent paper,i Chen and Jayaprakash studied a 
Kondo impurity with a pure power-law exchange cou- 
pling under conditions of strict particle-hole symmetry. 
An interesting feature of this work was a comparison be- 
tween the impurity contribution to the total ma gnetic 
susceptibility [the quantity Ximp defined in Eq. ( 5.11 ) 
above] and a local susceptibility xioc which directly 
probes the magnetic properties of the impurity. It was 
argued that the behavior of xi oc provides evidence for the 
existence of a finite critical coupling J c for r > \ — even 
though the thermodynamics provide no sign of such a 
critical coupling — and that the low-temperature physics 
is fundamentally different for J < J c and Jo > J c . In 
this subsection we re-examine this issue, and conclude 
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FIG. 17. Low-temperature limit of the local impurity sus- 
ceptibility Txioc vs the dimensionless Kondo coupling poJo 
in the exactly screened Kondo model, calculated for different 
values of the exponent r describing the power-law exchange. 
The data points were computed using a discretization param- 
eter A = 9. Solid lines are provided as a guide to the eye, 
and vertical dashed lines indicate the critical couplings for 
the cases r — 0.3 and r — 0.4. 
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FIG. 16. Schematic renormalization-group flow dia- 
grams for the exactly screened Kondo model, showing the 
Jo— Vo plane for Vi = and for fixed r: (a) < r < r*; 
(b) r* < r < r max ; (c) r > r max . Thin lines with arrows show 
the renormalization of effective couplings with decreasing 
temperature. Filled circles indicate RG fixed points. Broken 
lines represent flows out of the plane towards the asymmetric 
strong-coupling (ASC) fixed point, located at |Vi| = oo. 



that there is no such critical point for r > i, and that 
the local properties and total properties are perfectly con- 
sistent with one another. 

The local susceptibility xioc, denoted ((s z ;s z )) in 
Ref. H is a zero-frequency response function defined by 
the relation 



where s z is the z-component of the impurity spin; h is 
a magnetic field which acts only on the impurity and 
which enters the Kondo Hamiltonian through an addi- 
tional term ATix = g^Bhs z . 

In the standard case r = IX. xioc closely tracks Ximp 
as a function of temperature. C-2 Chen and Jayaprakash's 
results for r > can be summarized as follows: 

(i) For r < \ and J < J c (r), lim T ^ TXim P = 1/4 
while Txioc approaches a smaller, but still nonzero, value. 

(ii) For r < \ and J > J c {r), Tximp beads to the 
value r/8, whereas lim/r^o ^Xioc = 0. 

(iii) For r > | there exists a finite critical value, 
paJ c = 0(r), such that the behavior for Jo < J c is the 
same as in (i). 

(iv) For r > | and Jq > J c , Txj mp approaches 1/4 but 
liniT^o Txioc = 0. 

The authors of Ref. [| interpreted result (iv) as in- 
dicating that the impurity spin is locally quenched 
even though the total magnetic susceptibility shows no 
Kondo effect. This implies the existence of a third low- 
temperature regime, in addition to those governed by 
the local-moment fixed point [reached in cases (i) and 
(iii)] and the symmetric strong-coupling fixed point [the 
ground state for (ii)]. 

Since the total thermodynamic properties for r > j 
seem to indicate flow to the local-moment fixed point 
for all Jo, we have systematically examined the zero- 
temperature limit of Txioc- Within the numerical RG 
framework the right-hand expression in Eq. (6.10) can 
be evaluated for a small but finite value of h. (This is 
the same method as employed by Chen and Jayaprakash. 
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FIG. 18. Same as Fig. [17], restricted to r < | and plotted 
on a log-log scale. Data are shown for two different discretiza- 
tion parameters, A = 3 and 9. Straight lines represent fits 
to those A = 9 data points having an exchange coupling Jo 
within 3% of the critical value J c . 



and 0.7 offer no hint of critical behavior within a range of 
PoJo extending well beyond r. For large Jo, liniT^o ^Xioc 
falls off in roughly exponential fashion. It appears prob- 
able that for all Jo, Txioc heads to a nonzero value, de- 
scribing incomplete quenching of the impurity spin. 

These observations are consistent with our previ- 
ous conclusion,!!] based on the finite-size spectrum and 
the computed thermodynamic properties, that the low- 
temperature behavior of the particle-hole-symmetric 
Kondo model for r > \ is described quite straightfor- 
wardly by the weak-coupling (local-moment) limit. For 
r < i , the approach of Jo to its critical value from below 
is signaled by the vanishing of Txi oc , while the divergence 
of xi oc marks the approach from above. The behavior of 
Tx\oc makes it a candidate order parameter for describ- 
ing the J c critical point-. Further investigation of this 
possibility is under way.c£l 

3. Relationship between the Kondo and Anderson models 



It should be noted that for h ^ 0, the total sp in S is no 
longer a good quantum number.) Equation (3.1C) im- 
plies that a finite limiting value of Txioc reveals itself 
in a T _1 variation of (s z ). The field h must be chosen 
sufficiently small that any such T -1 regime can be de- 
tected before (s z ) saturates, as it must ultimately do 
since \s z \ < \. Our calculations were performed for 
10 -13 < gpsh/D < 10~ 7 . As a check, it was estab- 
lished that the results were insensitive to the precise field 
used. In order to reduce the computer time required for 
these runs, we used A = 9 with E c = 40. As mentioned 



in Section VC, it is possible to obtain accurate zero- 



temperature results even with such a large discretization 
parameter. We have verified that reducing A may pro- 
duce small shifts in the critical coupling J c , but does not 
change our essential conclusions (presented below). 

Figure [l7] plots the low-temperature limit of T\\oc 
versus the exchange coupling pqJq for four values of r. 
Consider first the data for r — 0.3 and 0.4, which in- 
dicate that limT_»o ^Xioc is finite for all Jo < J c , and 
vanishes for Jo > J c . Here, the critical coupling J c 
coincides with that deduced from the thermodynamic 
properties or the low-temperature many-body eigenspec- 
trum: p J c (r = 0.3) = 0.343 and po J c (0.4) = 0.491 (both 
to three significant figures). As shown in Fig. |lj, 
the form of the curves for Jo just below J c is well- 
described by a power law, limT_>o Txioc oc (J c — J) ly , 
with v[r = 0.3) = 0.70 ± 0.05 and u(0A) = 1.80 ± 0.08. 
(Figure [l8| also presents data for r — 0.4 computed us- 
ing a discretization parameter A = 3. The resulting 
exponent, v = 1.81 ± 0.05, is in close agreement with 
that obtained using A = 9.) On the strong-coupling 
side of the critical point (not shown) it is found that 
lim T ^oXiocai oc (J - J C )~ A , where A(r = 0.3) = 4.3 ±0.1 
and A(0.4) = 3.9 ±0.1. 

By contrast, the curves in Fig. [Tt] representing r = 0.6 



The mapping between the Anderson and Kondo mod- 
els via the Schrieffer- Wolff transformation [Eqs. ( 2.10 )] 
can be formally justified only for — e^, U 3> T and 
< poJo, | A) Vb| *C 1- In these limits, the physical prop- 
erties of the Kondo mode l must be equivalent to those 
presented in Section VI A provided that one equates the 



half-bandwidth D (Kondo) with Tp (Anderson), the tem- 
perature at which charge fluctuations freeze out. For the 
standard case r = 0, the properties of the Kondo model 
are universal for all Jo > and can be fully explored 
from within the Anderson model. The relation between 
the two models becomes nontrivial, however, in the pres- 
ence of a power-law scattering rate, which introduces a 
new energy scale J c into the problem. 



For small, positive values of r, Eqs. (2.10) can yield 
values of Jo greater than J c , as well as values less than 
J c . One therefore expects all behaviors of the Kondo 
model to be reproduced by the Anderson Hamiltonian. 
This is confirmed by Fig. [l^, which superimposes ther- 
modynamic properties computed for the particle-hole- 
symmetric Kondo and Anderson models with r = 0.2. 
The figure shows flows to weak coupling and to sym- 
metric strong coupling, as well as the critical behavior 
associated with the intermediate-coupling fixed point. 
Each Anderson curve — computed for one of the param- 
eter sets from Fig. pi and plotted versus T/\ed\ because 
Tp ~ \ed\ in these examples — is reproduced almost ex- 
actly by the Kondo model with a suitable choice of Jo- 

For r > 0.5, by contrast, it is impossible to attain val- 
ues of the exchange Jo > J c under the Schrieffer- Wolff 
transformation, and the strong-coupling behavior of the 
Kondo model (e.g., the particle-hole-asymmetric thermo- 
dynamic properties in Figs. 2 and 3 of Ref. 0) cannot be 
reproduced. This does not mean that strong coupling is 
unattainable within the Anderson model; only that the 
route to strong coupling via the local-moment regime 
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FIG. 19. Impurity susceptibility Tximp and entropy 
Simp/fe for r = 0.2, plotted vs either T/D (screened Kondo 
model with Vo = 0, solid symbols) or T/|eti| (symmetric An- 
derson model with To = 0.1D, open symbols). In four of the 
six cases shown, the entropy rapidly converges with decreas- 
ing temperature to the value ks In 2. 



is blocked. A direct crossover from valence fluctuation 
to asymmetric strong coupling (or, cquivalently, to the 
frozen-impurity regime) takes place for all r > when- 
ever the impurity level is placed above, or only slightly 
below, the Fermi energy. 

We now turn to a comparison between the 
intermediate-coupling fixed points of the Kondo and An- 
derson models. A number of pieces of evidence point to 
the complete equivalence of these fixed points: (1) The 
thermodynamic properties coincide to within the accu- 
racy of our calculations. Figu re fl9| illustrates this ex- 
plicitly for r — 0.2, while Figs. |10] and 14 show that for 

Similar agree- 



0.2, while Figs. 
r = 1, Tximp ~ 0.164 in both models, 
ment is found for other values of r, as well. (2) An 
extensive comparison of the energies and quantum num- 
bers of the many-body eigenstates indicates that the low- 
energy spectra at the Kondo and Anderson fixed points 
are identical. In each model, there is a range of expo- 
nents, r* < r < r max , over which there are two distinct 



intermediate-coupling fixed points (see Section [VI B l| ) . 
Within this range, the spectra at the particle-hole sym- 
metric fixed points of the two models match, as do the 
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FIG. 20. Position of the intermediate-coupling fixed point 
poJc (defined in the text) vs r for the Anderson model with 
a pure power-law scattering rate. Solid lines are provided 
as a guide to the eye. The dashed line shows the correspond- 
ing quantity computed for the particle-hole-symmetric Kondo 
model (data from Ref. Q). 



levels at the nonsymmctric fixed points. (3) It is shown in 
the next paragraphs that for small r, the critical coupling 
r c of the Anderson model corresponds to an effective ex- 
change coupling very close to the value of J c measured 
directly in the Kondo model. (We argue below that the 
positions of the two fixed points are not expected to co- 
incide for larger values of r.) 

The position of the two intermediate-coupling fixed 
points can be compared using the scaling theory dis- 
cussed in Section VI A3. Equations (EO)-(6.5) and (|6.9|) 



may be combined to convert the critical value of the bare 
scattering rate, r c plotted in Fig. |^, into the dimension- 
less exchange poJ c on entry to the local- moment regime. 
Figure ^0] shows that this transformation collapses the 
two data sets for U = -~2e<i onto a single curve. A dashed 
line shows the critical exchange computed directly within 
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the particle-hole-symmetric Kondo model.0e3 The criti- 
cal couplings of the two models are in close agreement 
for small r. The significant differences that develop for 
larger r, in a range where T c becomes of order (e^l, may be 
attributed to a breakdown of the Schrieffer- Wolff trans- 
formation as charge fluctuations become important. (In 
Rcf. |l^ a similar comparison is made between the effec- 
tive values of r c in the Anderson and Kondo models.) 

The infinite-?/ curves plotted in Fig. ^ lie close to each 
other, and to the data for a symmetric impurity, only for 
r < 0.25. For larger r, J C (T C ) turns sharply downward, a 
feature not seen in the particle-hole- asymmetric Kondo 
model (see Fig. 1 of Ref. 0). It is clear from Fig. || that 
the position of the critical point is much better described 
by a relation of the form T c oc r\ed\ than by poJ c <x r. 
Again, charge fluctuations on the impurity site account 
for these differences. The Kondo intermediate-coupling 
fixed point is always approached from within the local- 
moment regime. For large r, however, the T c fixed point 
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of the Anderson model is instead reached directly from 
the high-temperature regime (see Fig. [Io| for insta nce) . 
Indeed, comparison of with t* d _ given by Eq. (6.7) 
shows that the value of T c places the system outside 
the local-moment regime for all r > 0.25 in the case 
€d = —0.1D, and for all r > 0.15 in the case = — 0.01Z). 
We conclude, therefore, that in all situations where the 
mapping of the Anderson model onto the Kondo model 
is justified, the position of the intermediate fixed point 
determined in the two problems is in good agreement. 

In summary, the RG fixed points of the Kondo model 
appear to form a true subset of those of the Anderson 
model. For r > 0.5, though, certain paths between these 
fixed points that can be followed in the pure-spin problem 
cannot be realized once charge fluctuations are allowed. 
In this sense, the Kondo model with power-law scattering 
has an existence independent of the Anderson model. 



C. Underscreened Kondo model 

In this section we present results for the Kondo Hamil- 
tonian describing the interaction of a spin-one impurity 
with a single band of electrons. We focus on the differ- 
ences between this model and the conventional (exactly 
screened) problem addressed in the previous section. 

The weak-coupling limit of the s — 1 model shares 
many features with the conventional case. In particular, 
the fixed-point at Jo = is marginally u nstab le for r = 
but (according to the analysis in Section IV D ) it is stable 
for all r > 0. For small positive values of r, one can apply 
poor-man's scaling to demonstrate!!!] the existence of an 
intermediate-coupling fixed point at po J c ~ v. 

At antiferromagnetic strong coupling, the s = 1 model 
behaves very differently from the s — ^ problem. The 
larger impurity spin is "underscreened" and retains a net 
spin r = i . For r = this limit is known to be marginally 
stable^! but (as shown in Section IV E] ) for any r > the 
residual impurity spin destabilizes the symmetric strong- 
coupling fixed point. 

Given the instability of the fixed points at Jo = J c and 
Jo = oo, one might expect any Jo > J c to produce flow to 
asymmetric strong coupling. At particle-hole symmetry 
this option is ruled out, however, suggesting the existence 
of a stable fixed point at some exchange coupling that we 
shall denote J*, where J c < J* < oo. Our numerical RG 
calculations support this conjecture, at least for small r. 

We begin by considering the particle-hole symmetric 
case, Vq = 0. Figure ^l] shows the impurity susceptibil- 
ity computed for r = 0.2. For all values of the exchange 
coupling Jo , Tximp varies rather slowly with temperature 
and heads to a nonzero T = limit. However, closer in- 
spection reveals three qualitatively distinct cases. Tun- 
ing poJo to 0.254 yields a flat curve reminiscent of the 
(unstable) J e fixed point of the exactly screened model 
(see Fig. |l9|), whereas any smaller exchange produces a 
monotonic rise in T\imp towards the weak-coupling value 
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FIG. 21. Impurity susceptibility Tyi mp vs temperature T 
for the s = 1 Kondo model with pure power-law scattering 
specified by r = 0.2. 



of 2/3. For poJo > 0.254, Tximp falls towards a low- 
temperature limit of approximately 0.33, significantly 
higher than the strong-coupling value, (2+ri)/8 = 0.275. 
We interpret this as one piece of evidence for the J* fixed 
point alluded to in the previous paragraph. 

Behavior qualitatively similar to that shown in Fig. |2l| 
is found for all < r < r max , where 0.26 < r max < 0.27. 
For all r > r max , by contrast, the thermodynamics 
contain no signature of any intermediate-coupling fixed 
point. Instead, the system flows to weak coupling 
(Tximp = 2/3) for all values of the bare exchange J . 

Some of the systematic trends with increasing r are 
shown in Figs. [22]-^4[ Figure |H plots (filled symbols) 
the position po J c of the unstable intermediate-coupling 
fixed point as a function of r. For small r, po Jc ~ r, as 
expected from poor-man's scaling, but the curve turns 
upward and then abruptly terminates at r = r max . 

The critical coupling J* corresponding to the second 
intermediate-coupling fixed point cannot be determined 
directly because the low-temperature behavior does not 
undergo any qualitative change as Jo passes through J*. 
However, one can estimate the value of J* by examining 
the low-energy many-body spectrum. The eigenvalue of 
the first excited state at each of the four fixed points of 
the symmetric underscreened problem (Jo = 0, J c , J* 
and oo) is plotted versus r in Fig. ^3[ For small r, the 
deviation of the J c levels from the weak-coupling ener- 
gies is linear in r, consistent with the relation po J c ss r. 
It appears that J* is infinite at r — (at which point 
the extrapolated levels for the stable intermediate fixed 
point coincide with the strong-coupling values), and falls 
progressively as r increases. The energies at the two in- 
termediate fixed points are projected to cross at a value 
of r between 0.26 and 0.27 — precisely the range in which 
these fixed points disappear. 

Figure shows the impurity susceptibility and en- 
tropy at the two intermediate-coupling fixed points. For 
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FIG. 22. Critical couplings vs r for the s = 1 Kondo model 
with pure power-law scattering. Solid lines are provided as 
a guide to the eye. The particle-hole-symmetric fixed-point 
coupling J c (filled symbols) is plotted both to scale and mag- 
nified x5. For nonzero potential scattering (open symbols) 
there can be zero, one, or two critical couplings (J' < J"), 
depending on the value of r. 



each property, the J c and J* curves deviate with increas- 
ing r from the weak- and strong-coupling limits, respec- 
tively. Just as for the eigenenergies, the two sets of curves 
appear to cross in the range 0.26 < r < 0.27. Both 
^Ximp and vary linearly for small r. Empirically, 
the susceptibilities for r < 0.15 are well-described by the 
formulae 

TXim P ( Jc) = | - 0.7r, T Ximp (r) = \ + 0Ar, (6.11) 

while the entropies for r < 0.2 have an excellent fit to 
the weak- and strong-coupling expressions: 



S imp (J c ) = In 3, S imp (J*) = In 2 + rln4. 



.12) 



We have been unable to reproduce these properties by 
constructing suitable phase shifts for noninteracting elec- 
trons. We speculate that these intermediate-coupling 
fixed points are truly non-Fermi-liquid in character. 

Let us summarize the situation at particle-hole sym- 
metry. For < r < r max there are two intermediate- 
coupling fixed points — one stable (J*), the other un- 
stable [J c )- These fixed points merge at r = r max , above 
which value they both disappear, leaving weak coupling 
as the only stable fixed point. (We remind the reader that 
in the exactly screened Kondo model, the disappearance 
of the J c fixed point can be tied directly to the value 
r ma i = i at which the symmetric strong-coupling fixed 
point becomes unstable. We have no such argument to 
fix the precise value of r max in the underscreened model.) 

The introduction of potential scattering considerably 
modifies the picture presented above. As shown in Sec- 



FIG. 23. Lowest eigenvalues E* of the discretized a — 1 
Kondo model (A = 3) plotted vs r. Data are shown for the 
four fixed points of the particle-hole-symmetric problem, both 
for TV even (open symbols) and for N odd (filled symbols). 
Solid lines are provided as a guide to the eye. The J c and 
J* curves are extrapolated at r = to the weak-coupling and 
strong-coupling values, respectively. 



destabilizes the J* fixed point — found in the symmet- 
ric problem for all < r < r max — towards asymmetric 
strong coupling. The effect on the J c fixed point is more 
subtle, as will be explained in the paragraphs that follow. 
There are several parallels with (but also clear differences 
from) the behavior of an exactl y screened Kondo impu- 



rity described in Section VI B 1 



tion IV D, this perturbation is irrelevant in the weak- 
coupling regime. However, our numerics indicate that it 



We first consider values of r less than r max . For each 
r in this range and for any Vq ^ one can find a critical 
coupling J'(r,Vo) such that any Jo < J' yields weak- 
coupling behavior while any Jo > J' drives the system to 
asymmetric strong coupling. (The latter case is exempli- 
fied by the curve for r — 0.2, p J — 0.45, and p Vo = 0.3 
in Fig. ^l], which shows T^imp heading towards its asym- 
metric strong-coupling value of 1/4.) 

The behavior for Jo precisely equal to J'(r, Vq) de- 
pends on the value of r. For < r < r* sa 0.245, particle- 
hole asymmetry is irrelevant on the separatrix, and the 
system approaches the fixed point located at Vq = 0, 
Jq = J c = J'(r, 0). The RG flows for this case are 
sketched in Fig. |25|(a). For r* < r < r max , by con- 
trast, the flow along the line Jo = J'{r, Vq) is towards 
a new intermediate-coupling fixed point at Vq = V c , 
Jo = J'Jt) = J'( r >Vc), as shown in Fig. ^5|(b). (As was 
the case for a screened impurity spin, there is actually 
a pair of J' c fixed points at Vq = ±V C . Throughout this 
section, these two fixed points — which are related by 
particle-hole interchange — are treated as one, and all 
properties discussed will be assumed to depend only on 
the absolute value of Vo.) 

For r > r max , there is a range of potential scatter- 
ings | Vo| < V c over which the low-temperature physics 
is governed by the weak-coupling fixed point, whatever 
the bare exchange coupling J , i.e., no critical exchange 
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(a) 



FIG. 24. Impurity susceptibility Tximp and entropy S mp 
at the two intermediate-coupling fixed points of the s = 1 
Kondo model, plotted vs r. Solid lines are provided as a 
guide to the eye. The J c and J* curves are extrapola ted a t 
r = to the values for w eak coupling [given by Eq. (5.24)] 
and strong coupling [Eq. (5.34)], respectively. 



J'(r, Vq) can be found. For |Vo| > V Cl by contrast, there 
are two critical couplings. The system flows to weak cou- 
pling both for J < J'(r,Vo) and for J > J"(r,Va), 
while asymmetric strong coupling is reached for cou- 
plings which fall between these critical values. If Jo is 
tuned precisely to J' or to J", the system flows to an 
intermediate-coupling fixed point which we take to be lo- 
cated at V = V c , J = J' c (r) = J'(r, V c ) = J"{r, V c ). The 
RG flows that we deduce for this range of r are sketched 
in Fig. H(c). 



It was shown in Section [V E that there are two equally 



relevant operators in the vicinity of symmetric strong 
coupling (Jo = oo, Vq = 0). Figure |25| illustrates the com- 
petition between Oj 1 , which drives the system towards 
weak coupling, and Ov 1 , which causes flow towards asym- 
metric strong coupling. For r < r max , the J c and J* 
fixed points block flow along the axis Vq — 0, allow- 
ing potential scattering to dominate the low-temperature 
behavior. For r > r max , by contrast, the sole surviving 
intermediate-coupling fixed point at (J' c , V c ) stifles the 
growth of particle-hole asymmetry and instead steers the 
system to weak coupling. 

We now present some of the numerical evidence in 
support of the picture laid out above. Figure ^|| plots 
J' for fixed p Vo and for < r < 0.7. For r > r max , 
the second critical coupling J" is also plotted. One sees 
that J" is typically very large (greater than the band- 
width), implying that the upper bound on the asymmet- 
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FIG. 25. Schematic renormalization-group flow diagrams 
for the s = 1 Kondo model, showing the Jo—Vo plane 
for Vx = and fixed r: (a) < r < r*; (b) r* < r < r mw ; 
(c) r > r max . See Fig. |l(] for an explanation of the symbols. 



ric strong-coupling regime is unlikely to be accessible in 
practice. We identify the rightmost point on each crit- 
ical curve (the meeting of the J' and J" curves) with 
the J' c fixed point. Note that this particular plot yields 
r and poJ' c as functions of poV c . However, by invert- 
ing the procedure to make r the independent variable, 
one can deduce that both J' c and V c are increasing func- 
tions of r, at least over the parameter range shown. One 
can deduce, for instance, that poV c (r = 0.38) w 0.1 and 
PoV c (r = 0.68) w 0.5. These rather large values of V c sug- 
gest that for r 0.5, flow to asymmetric strong coupling 
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FIG. 26. Energy splittings AE* of the lowest pair of 
charge-conjugate eigenstates at the particle-hole-asymmetric 
intermediate-coupling fixed point of the s = 1 Kondo model, 
plotted vs the reduced exponent r/r* — 1, for A = 3 and 
r* = 0.245. Straight lines show fits to the leftmost four data 
points for odd- and even-numbered iterations of the numeri- 
cal RG method. 



(i.e., Kondo-screening of the impurity) can be achieved 
only under conditions of strong particle-hole asymmetry. 

Figure ^6] shows the energy splitting AE* of the low- 
est pair of charge-conjugate states at the J' c fixed point 
for 0.25 < r < 0.6. For r < 0.3, this splitting is well- 
fit by AE* cx (r - r*f\ with r* = 0.245 ± 0.002 and 
v* =0.44 ±0.15. For r < r* the fixed-point spectrum 
is observed always to be particle-hole-symmetric. 

The impurity susceptibility and entropy at the vari- 
ous intermediate-coupling fixed points are compared in 
Fig. ^4| Over the rather narrow range r* < r < r max in 
which the J c and J' c fixed points coexist, their properties 
are seen to diverge steadily. For larger r, both Txim P (J' c ) 
and Si mp (J' c ) rise monotonically. 

Finally, we note that all properties of the J' c fixed point 
shown in Figs. |2^, ^J, and ^6] appear to vary smoothly 
as r passes through r max , even though the fixed point is 
defined in a different manner for r < r max and r > r max . 
This serves as an indication that the fixed-point couplings 
( J' c , V c ) evolve continuously across the border between 
the regimes shown in Figs. 123(b) and (c). 



D. Overscreened Kondo Model 

In this section we present results for the Kondo Hamil- 
tonian describing the interaction of a spin-one-half im- 
purity with two degen erate bands or channels of elec- 
trons [Eqs. ( 2.12 ) and ( 2.13] )]. The weak-coupling prop- 
erties of the two-channel problem are very similar to 
those presented in the previous two subsections — includ- 
ing the existence for small r of an unstable fixed point 
at pqJ c ~ r. We focus on the intermediate-to-strong- 
coupling regime, where the three models differ markedly. 



At strong coupling, the impurity is "overscreened" by 
the two conduction bands and retains a net spin of. one- 
half. For r — 0, this limit is marginally unstable,Eil giv- 
ing rise to a stable, intermediate-coupling fixed point at 
PqJ* — 0(1). For r > 0, the symmetric strong-coupling 
fixed point is outright unstable due to the residual impu- 
rity degree of freedom, so one might again expect flow to 
some J* > J c , just as in the underscreened model. 

The existence of two intermediate-coupling fixed points 
was predicted previously using an extension of With- 
off and Fradkin's poor-mark's scaling analysisa to the 
A^-channel Kondo problem.!! The dimensionless Kondo 
coupling PqJ was found to rescale from its bare value 
PqJq according to the equation 



d(pgJ) 
dlnT 



= rp a J - (po-J) 2 + c{p jf 



(6.13) 



The coefficient c is a complicated function of r and N c , 
which reduces to c = N c / 2 in the limits r <C 1 and 
N c ^> 1. Then Eq. ( 6.13 ) has fixed points satisfying 
d{p J)/dlnT = at p J = 0, oo, and (1 ± y/1 - 2N c r). 
For small r, the intermediate-coupling fixed points are 
located at pqJ c ~ r (unstable) and po J* w 2/N c — r 
(stable). However, these two fixed points merge at 
r = l/2N c , p J=l/N c . For r > 1/2N C there is no 
intermediate-coupling fixed point and the RG trajecto- 
ries flow directly from strong coupling to weak coupling. 

Although the two-channel case does not strictly satisfy 
the condition N c ^> 1, the predictions of scaling theory 
are well borne out by numerical RG calculations. More- 
over, the intermediate fixed points turn out to survive 
the inclusion of potential scattering (which was not taken 
into account in Ref. [|). 

The addition of a second conduction band greatly in- 
creases the size of the basis of the discretized version 
of the Kondo model, and hence the computer time re- 
quired for a solution of the problem. We have found 
it impractical to compute thermodynamic properties us- 
ing a discretization parameter A = 3 while keeping all 
states up to an energy cutoff E c > 25, as was done in 
the single-channel problems. This prevents reliable de- 
termination of the temperature-dependence of the ther- 
modynamic properties. However, experience indicates 
that fixed-point properties can be computed accurately 
for values of A as large as 10. Below we present many- 
body eigenstates computed for A = 3, but thermody- 
namic properties obtained using A = 9 and E c = 25. 

Figure [27] plots the position poJ c of the unstable 
intermediate-coupling fixed point as a function of r. For 
small r, poJ c ~ r, as expected from poor-man's scal- 
ing. At particle-hole symmetry, the curve turns up- 
ward and then abruptly terminates at r = r max where 
0.23 < r max < 0.24. In this respect, the underscreened 
and overscreened problems are very similar. 

Just as for the underscreened problem, one can obtain 
indirect evidence for the value of J* by examining the 
low-lying many-body spectrum. The energy of the first 
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FIG. 27. Critical couplings for the two-channel Kondo 
model with pure power-law scattering. Solid lines are pro- 
vided as a guide to the eye. The data for poVo = and 
0.1 are plotted both to scale and magnified x4. For nonzero 
potential scattering there can be zero, one, or two critical 
couplings (J' < J"), depending on the values of r and Vo. 

excited state at each of the four fixed points of the sym- 
metric overscreened problem (Jo = 0, J c , J* and oo) is 
plotted versus r in Fig. The levels at the J c fixed 
point progressively diverge from the weak-coupling ener- 
gies as r is increased from zero, consistent with the re- 
lation poJ c w r reported in the previous paragraph. The 
smooth evolution of the levels at the stable intermediate 
fixed point suggests that J* decreases continuously with 
increasing r from its value for a constant scattering rate, 
PoJ*{r = 0) w 1. The energies at the two intermediate 
fixed points can be extrapolated to cross in the range 
0.23 < r < 0.24 where both fixed points disappear. 

Figure ^9] shows the impurity contributions to the 
susceptibility and the entropy at the two intermediate- 
coupling fixed points. For each property, the J c and J* 
curves deviate with increasing r from their weak-coupling 
and r — non-Fermi-liquid values, respectively, and the 
two curves can be extrapolated to cross in the range 
0.23 < r < 0.24. Both Txi mp and S[ mp vary linearly 
with r for r < 0.15. Empirically, the susceptibilities are 
well-described by the formulae 

Tximp(Jc) = \ - \, T Ximp (r) = T -, (6.14) 

while the entropies fit 

Simple) =ln2, S imp (J*) = im2 + rln4. (6.15) 

Two factors greatly impede the study of the effects of 
potential scattering in the overscreened Kondo model: 
(1) Away from part icle- hole symmetry, the total axial 
charge (see Section |IIF| ) is no longer a good quantum 
number. This change roughly doubles the size of the basis 
at each iteration of the numerical RG procedure. Even 




r 



FIG. 28. Lowest eigenvalues E* of the discretized 
two-channel Kondo model (A = 3) plotted vs the exponent r 
describing the power-law exchange. Data are shown for the 
four fixed points of the particle-hole-symmetric problem, both 
for N even (open symbols) and for N odd (filled symbols). 
Solid lines are provided as a guide to the eye. The J c curves 
are extrapolated at r = to the weak-coupling values. 

working with a discretization parameter as large as A = 
9, we have found it feasible to retain only those many- 
body eigenstates with scaled eigenvalues E* < E c m 15 
(compared to E c — 25 for the particle-hole-symmetric 
problem). (2) The instability of the two-channel Kondo 
model with respect to channel asymmetry^ is found to 
rise markedly with increasing r and \Vq\. Over much of 
the parameter space, unavoidable numerical asymmetry 
at the level of the machine precision grows to of order 
unity before the many-body energy levels get close to the 
zero-temperature fixed point of the channel-symmetric 
problem. In light of these obstacles, we focus our remarks 
on the qualitative features of the RG flow diagrams for 
r < 0.3 and \p V \ < 0.5. 

For all < r < r max , there appears to be a critical cou- 
pling J'(r, Vq) for any potential-scattering strength Vq. 
Figure |27] plots J'(r) for two fixed values of Vo- For 
Jo < J , the system flows to weak coupling, while for 
Jo = J' it reaches the J c fixed point of the particlc-holc- 
symmetric problem. For Jo > J', the flow is to a gener- 
alization of the J* fixed point found for Vo = 0. Specif- 
ically, the energy levels are obtained from those of the 
J* fixed point by splitting each pair of charge-conjugate 
states. At fixed Jo, this splitting grows with increasing 
Vb; with Vo held fixed and Jo starting at J', the split- 
ting initially grows from zero as Jo increases, then passes 
through a maximum, and eventually falls back towards 
zero as Jo — + oo. From this behavior, we deduce that the 
RG flows have the form shown in Fig. |30"|(a). To within 
the accuracy that we can achieve (around 2%), the im- 
purity entropy is the same everywhere along the line of 
fixed points, but T^imp falls as one moves away from the 
symmetric fixed point (see the inset to Fig. G9). 

The RG flows for r > r max [Fig. |30|(b)1 are qualita- 
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FIG. 29. Impurity susceptibility Tximp and entropy 
Simp at the two intermediate-coupling fixed points of the 
two-channel Kondo model, plotted vs r. Solid lines are pro- 
vided as a guide to the eye. The J c curves are extrap olate d 
at r — to the values for weak coupli ng, g iven by Eq. (5.23). 
The strong-coupling properties [Eq. (5.35)] are also plotted. 
Inset: Impurity susceptibility (in the same units as the main 
figure) at the J* fixed point vs potential scattering poVo, for 
r = 0.2 and po Jo = 0.6. 




FIG. 30. Schematic renormalization-group flow diagrams 
for the two-channel Kondo model, showing the Jo-Vb plane 
for Vi = and fixed r: (a) < r < r max ; (b) r > r max . See 
Fig. for an explanation of the symbols. 



tively similar to those of the underscreened s = 1 Kondo 
model. There exists a single intermediate-coupling fixed 
point at J = J 'J?), V = V c (r). For all \V \ < V c the 
system flows to weak coupling, whatever the value of the 
bare exchange coupling J . For \Vq\ > V c , by contrast, 
the flow is to weak coupling for Jo < J'(r,Vo) and for 
Jo > J"(r, Vq); otherwise the system flows to asymmet- 
ric strong coupling. Figure |27] plots J'{r) and J"(r) for 
PoVo — 0.5. The absence of any critical coupling for 
PoVo = 0.1 beyond r = r max indicates that poV c (r) > 0.1 
for all r max < r < 0.4. This observation and the discon- 
tinuity in the slope of the J'(r) curve for p Vo = 0.5 both 
suggest that the properties of the J' c fixed point do not 
vary smoothly as r passes through r max . 

We note that, unlike the exactly screened and under- 
screened models, the two-channel Kondo model does not 
seem to exhibit any range of exponents r* < r < r max 
within which particle-hole-symmetric and asymmetric 
versions of the J c fixed point coexist. Due to the nu- 
merical difficulties mentioned above, this possibility can- 
not be completely ruled out, but any range of dual fixed 
points is certainly very narrow (r max — r* < 0.02). 



E. Departures from a pure power-law scattering rate 

Now we consider various changes to the form of the 
power-law scattering rate defined in Eq. (1.3). We focus 
on three features which are likely to be present in real 
materials: (i) removal of the symmetry r(— e) = T(e); 
(ii) restriction of the power-law variation in T(e) to a re- 
gion of halfwidth A -C D; and (iii) the existence of a 
small but nonvanishing scattering rate at the Fermi en- 
ergy, r(0) ^ 0. The effects of these modifications can be 
predicted qualitatively using poor-man's scaling, and can 
be investigated in detail via numerical RG calculations. 
In the latter approach, each change in the for m of r(e) 
simply alters the values of F defined in Eq. ( |2.5| ) and 
the tight-binding coefficie nts e„ and t n entering the dis- 
crctized Hamiltonian, Eq. (2.32); otherwise the numerical 
treatment remains the same as for a pure power-law scat- 
tering rate. It turns out that the first modification above 
is relatively inconsequential, whereas the second can sig- 
nificantly increase the likelihood of observing the Kondo 
effect for values of r > |, and the third can produce even 
more fund ame ntal departures from the results obtained 
using Eq. ( Q . 
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1. Particle-hole- asymmetric scattering rate 



It was pointed out in Section III that if the scattering 



rate satisfies r(— e) = T(e) for all e, then the param- 
eters e n entering Eq. (2.32) are identically zero. This 
symmetry is unlikely to be exactly preserved in any real 
system. We have studied the effect of various symmetry- 
breaking per turb ations, the simplest being the modifica- 
tion of Eq. (1.2) so that the conduction band extends 
over energies —(1 + p)D < e < (1 — p)D. Setting /i ^ 
changes not only the hopping coefficients t n entering the 
discretized conduction band, but also leads to nonzero 
diagonal coefficients e„. This invalidates the asymp- 
totic expressions for the conduction-band eigenvalues and 
eigenvectors presented in Section [II, and vitiates the 



fixed-point stability analysis of Section IV. However, nu 



merical RG results indicate that the physical behaviors 
of the Anderson model remain essentially the same. The 
band asymmetry prevents the appearance of the sym- 
metric strong-coupling fixed point, but the other regimes 
discussed in Section exist and moreover exhibit the 
same power-laws in their thermodynamic properties. The 
effect of band asymmetry on the intermediate-coupling 
fixed points of the three variants of the Kondo model is 
essentially equivalent to that of potential scattering on 
the same problem with a symmetric band. 



2. Restricted power-law scattering rate 

In real gapless systems, the power-law variation of the 
scattering rate or of the Kondo exchange is unlikely to 
extend over the entire band in the manner assumed in 




Eq. (1.3). A more realistic form for the scattering rate 
rolls over to a roughly constant value outside a pseudogap 
of halfwidth A. This can be approximated by writing 



(6.16) 



The effects of restricting the power-law scattering 
regime can be predicted using poor-man's scaling. Con- 
sider a localized level described by the nondegenerate An- 
derson model. At temperatures T> A, the impurity is 
insensitive to the presence of the pseudogap, and one ex- 
pects the standard physics exhibited in a metallic host. 
As the temperature is lowered, the effective position of 
the impurity level scales upward according to Eq. (3.2), 
while To remains essentially constant. Two qualitatively 
different situations can arise. In the first, the impurity 
remains in the valence-fluctuation regime (see Fig. |lj) 
all the way down to temperatures T« A. In this case, 
once the temperature falls much below A, the system 
will behave very much like an Anderson impurity with 
a pure power-law scattering rate, the role of the half- 
bandwidth D being taken by A and with e</ replaced 



by £d(A) [given by Eq. (6^)]. The qualitative effects 
of the p seudo gap shoul d therefore be those reported in 
Sections VI A and VI B , although the magnitude of these 
effects will decrease as the pseudogap narrows. For exam- 
ple, on any subsequent entry to the local-moment regime, 
the effective Kondo exchange J will be reduced relative 
to the case r = by a factor of at least |e ( ;(A)/A| r ' (com- 
pared to a reduction of \ed/D\ r for A = D). 

Should there exist a solution to Eq. ( |6 . 5| ) such that 
Tp > A, then real charge fluctuations on the impurity 
site will be frozen out before the power-law density of 
states makes its presence felt. Perhaps the most inter- 
esting situation arises when the system enters the local- 
moment regime, in which case one can Schrieffer- Wolff- 
transform to the Kondo description of the problem. Over 
the temperature range T F > T > A, the conduction 
band will begin to screen out the impurity moment; in 
the scaling picture, the effective value of the Kondo cou- 
pling will renormalize upwards according to Eq. ( 6.13j ) 
with r = 0. At temperatures T <C A, the impurity 
maps onto a model with pure power-law exchange having 
an effective half-bandwidth A and an exchange coupling 
p J « [ln(A/T° )]-!. Here T% = Dexp[-l/(p J )}, the 
Kondo temperature for a system having a constant den- 
sity of states po, is assumed to be smaller than A. (If 
T^r > A, then the impurity will already have entered the 
strong-coupling regime before the pseudogap comes into 
play, in which event perturbative scaling can provide no 
insight into the behavior for T < A.) 

The scenario of the previous paragraph implies a sig- 
nificant enlargement of the region of parameter space 
within which a fully developed Kondo effect can take 
place. Consider, for instance, the case U — oo, with 
some fixed value t of the hybridization matrix ele- 
ment. Working to lowest order, let us neglect the many- 
body renormalization of the impurity energy (which 
beco mes in creasingly weak as r increases; see Sec- 
tion VI A 3). Then the exchange pqJq entering the ef- 
fective low-temperature Kondo problem is approximately 



PqJ° = t /(D\ed\) for a constant scattering rate 



for 



0), 
D), 



pure power-law scattering (A 

PoJ° m (l e d|/A)] for restricted power- 
The last value is enhanced over that ob 



PoJ°\e d /D\ r 

and po J° [l — 
law scattering, 
taincd with a constant scattering rate, and for A suffi- 
ciently small, Jo will exceed the threshold J c (r). This is 
true even for r > 0.5, a range in which no Kondo effect 
can be observed in cases of pure power-law scattering. 

The predictions made in the preceding paragraphs can 
be tested against numerical RG results. Restriction of 
the power-law scattering to a region lej < A alters the 
hopping coefficients t n entering Eq. ( |2.32 ). The values 
for small n [s uch that the characteristic temperature T n 



given by Eq. ( [4.3|) greatly exceeds A] become essentially 
identical to the corresponding values for the case r = 0, 
while the t n ' s for large n (such that T n <§C A) are still 
given by Eq. (3.3). Since it is the large- n coefficients that 
determine the low-temperature behavior of the system, 
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FIG. 31. Impurity susceptibility Tximp and entropy Si mp 
vs temperature T for the infmite-(7 Anderson mo del w ith a 
restricted linear scattering rate described by Eq. (3.16) with 
r = 1, A = 10" 3 D and T = 0.050D (symbols, solid lines). 
Data are also shown for the same impurity parameters (t and 
£d) but a constant scattering rate (dashed lines). 



to produce an upturn in both T^imp and S- lmp , signaling 
a weakening in the effective coupling between the impu- 
rity level and the conduction band. For e<2 > — 0.16D, 
this weakening is reversed as the temperature decreases 
further, and the system eventually flows to asymmet- 
ric strong coupling. In the case td = — Q.16D, where 
the impurity clearly enters the local-moment regime at a 
temperature well above A, the downturn in each prop- 
erty indicates the existence of a Kondo effect where none 
was found in the pure power-law case (see Fig. |l(]). For 
ed < —0.17D, the rise in Txi mp towards the free local- 
moment value indicates that when the system encoun- 
ters the power-law scattering regime, the effective ex- 
change coupling is smaller than the critical value needed 
to bring about Kondo-screening of the impurity. Some- 
where in the range —0.17 < ed/D < —0.16 there pre- 
sumably exists a critical impurity energy that places the 
system precisely at the intermediate-coupling fixed point, 
where T X i mp » 1/6 (see Fig. @). 

We have performed similar calculations for parame- 
ters chosen so that real charge fluctuations remain possi- 
ble when the temperature becomes comparable with the 
width of the pseudogap (i.e., for a larger A or a smaller 
|e<i| and To than used in Fig. |3l]). As surmised in the scal- 
ing discussion above, the physical properties depart less 
dramatically from those for a pure power-law scattering 
rate than in the case shown in Fig. [H]. 

We have also investigated the effect of restricted power- 
law exchange within the Kondo models. Results for the 
exactly screened case appear in Ref. [t]. In this and the 
other variants of the model, the low-temperature behav- 
ior can be accounted for by assuming a pure power-law 
exchange with the bare coupling Jo replaced by a (larger) 
effective value that compensates for the elimination of 
conduction band states at energy scales |e| > A. 



the analysis of the stable fixed points of the Anderson 
model in Section |^ remains applicable. 

Figure |3l] presents numerical RG data for the Anderson 
model with r = 1, A/D = 10~ 3 , and T /D = 0.050. [The 
value of To is chosen so that if the energy-dependence 
of r(e) were to derive solely from the density of states, 
then the (energy- independent) hybridization t would be 
identical to that used in Fig. nO.] The numerical results 
bear out quite well the predictions of poor-man's scaling. 
Over the range T > A, the impurity susceptibility and 
entropy are very close to those obtained using the same 
t and td but a constant scattering rate (dashed lines). For 
6d < — 0.16D, an initial increase in T\i mp signaling entry 
to the local-moment regime is followed by a downturn 
as the conduction band begins to screen the impurity 
moment; whereas for ed > —0.1D, Tx; mp falls monoton- 
ically in the manner characteristic of the mixed-valence 
and empty- impurity limits. 

The pseudogap in the scattering rate begins to make 
its presence felt about a decade in temperature above A. 
For all the values of Cd shown in Fig. Bl], the initial effect is 



3. Finite Fermi- energy scattering rate 

In the context of d-w&ve superconductivity, it has been 
predicted] that the pair-breaking effect of any finite con- 
centration of magnetic impurities will feed back to pro- 
duce a small but nonzero quasiparticle density of states 
at zero energy. This may be taken into account by in- 
troducing a lower cutoff A' on the power-law scattering 
rate (in addition to the upper cutoff A considered above) , 
such that r(e) « r (A'/A) r for |e| < A'. 

Based on a heuristic scaling approach, similar to that 
employed in the previous subsection, it is apparent that 
the system can exhibit a wide variety of behaviors, de- 
pending on the relative magnitudes of the energy scales 
in the problem. The lower cutoff on the power-law scat- 
tering rate cannot significantly affect the physics until 
the temperature drops to become of order A'. At this 
point, the clearest departure from the case A' = arises 
if the impurity is already in, or subsequently enters, its 
local-moment regime. In such cases the existence of a 
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nonzero Fermi-surface scattering rate ensures the even- 
tual appearance of a standard, metallic Kondo effect. 
However, for r > 1 at least, most plausible values of 
the various parameters in the model lead to screening of 
the impurity only at temperatures that are so low as to 
be physically irrelevant. 

VII. SUMMARY 

In this work, we have studied four models in which 
the interaction between a magnetic impurity level and an 
otherwise-uncorrelated host fcrmion system is described 
by an energy-dependent scattering rate that vanishes in 
power-law fashion at the Fermi level. Our principal re- 
sults are summarized in the following paragraphs. 

In all four models, the effect of the pseudogap in the 
low-energy scattering rate is to create a low-temperature 
regime which has no counterpart in metals (where the 
scattering rate is fiat in the vicinity of the Fermi level). 
This stable weak-coupling limit, first identified!] in the 
exactly screened Kondo model, is characterized by a free 
impurity local moment which retains its Curie suscepti- 
bility down to zero temperature. 

Depending on the value of r, the weak-coupling region 
of parameter space is bounded in each model by either 
one or two unstable intermediate-coupling fixed points. 
A particle-hole-symmetric fixed point of this type exists 
for all < r < r max and is the only such fixed point for 
< r < r* , over which range particle- hole asymmetry is 
marginally irrelevant. For r > r* a second fixed point 
with different properties is reached whenever particle- 
hole asymmetry is present, and for r > r max it is the only 
fixed point of this type. The values of r max and r* are 
model-dependent, but in all the four cases that we have 
studied < r* < r max < |. 

Another common feature of these models is the exis- 
tence of two distinct strong-coupling regimes. The first, 
symmetric strong coupling, is the natural generalization 
to r > of the low-temperature limit of the metallic 
Kondo problem. In the Anderson and screened Kondo 
models, the impurity degree of freedom is completely 
quenched, as evidenced by the vanishing of the quan- 
tity Txioc and the sole effect of the magnetic level is 
to impose a phase shift on electrons at the Fermi en- 
ergy. This phase shift results in anomalous thermody- 
namic properties, namely those that would arise as if a 
fraction min(r, 1) of an electron were to decouple from 
each conduction band. For exponents r < 1, the phase 
shift should also result in a zero-temperature peak in the 
electrical resistivity. The magnitude of this peak should 
decrease with increasing r, and for r > 1 it should disap- 
pear altogether. 

Except under artificial conditions of strict particle- 
hole symmetry, the symmetric strong-coupling limit is 
unstable. Generally, the system is driven to an asym- 
metric strong-coupling fixed point at which a Fermi-level 



phase shift of ±7r implies a vanishing impurity contribu- 
tion to thermodynamic and transport properties. Un- 
derscreened and overscreened impurities, however, admit 
two other possibilities: For r < 1/4, one can obtain flow 
to (or, for the s = 1 model, near) a non- Fermi-liquid fixed 
point, corresponding to an effective coupling J* which 
decreases with increasing r. For r > 1 /4, the generic be- 
havior is flow to weak coupling, in which case the absence 
of interesting many-body effects invites comparison with 
the ferromagnetic Kondo problem in metals. 

In addition to characterizing the location, stability, 
and thermodynamics of the various fixed points men- 
tioned above, we have also studied the possible crossovers 
between regimes dominated by these fixed points. We 
have quantified observations made previously based on a 
perturbative scaling analysis of the Anderson model that 
while local-moment formation is assisted by the presence 
of a pseudogap in the low-energy scattering rate, there 
is a strong suppression of the effective exchange coupling 
between any such moment and the delocalized electrons. 
As a result, it becomes progressively harder as r increases 
from zero — and essentially impossible for a pure power- 
law scattering rate with r > 0.5 — to find any choice 
of model parameters which leads to Kondo-screening of 
the impurity moment. In more realistic situations where 
the power-law scattering is confined to a narrow range of 
energies around the Fermi level, the suppression is less 
complete. Even in such cases, though, there is a signif- 
icant region of parameter space in which an unscreened 
local-moment survives down to zero temperature. 
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APPENDIX A: RESONANT LEVEL MODEL 
WITH A POWER-LAW MIXING RATE 

In this appendix, we analyze a noninteracting resonant 
level model for the mixing between an impurity level and 
a spinless conduction band. This model, described by 
the Hamiltonian 

H = e k 4ck + Qrf t rf + ^y=(c[<i + H.c), (Al) 

repre sent s the limit U = of the Anderson model 
[Eq. ( |2.l| )1, in which up and down spins decouple from 
one another, and the spin index can therefore be dropped. 
The mixing rate, 

r^Tr^i^e-ek), (A2) 

k 
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is as sum ed to have the pure power-law form given in 
Eq. (L3). Since the Hamiltonian is quadratic, various 
properties of the model can be calculated analytically. 



a. Impurity self-energy 

The impurity self-energy in the model described by 
Eq. flO) is 



EM = E 



The real and imaginary parts of this quantity are 



ReE(w) = P ^ 

k 

ImE(cj) = -r(w) 



l^kl 2 /^ _ 1 p 

U) — 6k 7T 



r(e) 



(A3) 

(A4) 
(A5) 



Substituting Eq. (1.3) into Eq. (A4), one obtains 



In 



D+LO 



D-uj 



0, 



ReS(w) = < 



-r(w)tan — sgn(w), 



r( w ) 



In 



£ 2 -cj 2 



2r w 



(r-l)7r £>' 



< r < 1, 
sgn(w), r = 1, 
r > 1. 

(A6) 



The expressions for r — and r = 1 are exact, but 
those for other values of r are approximations which are 
valid only for |w| < u) c (r). Here, u) c (r) is a cutoff energy 
scale which approaches zero as r approaches 1 from above 
or below. In all cases, Re EM passes through zero at 
uj = 0. However, ReE(w) has the same sign as u> for 
r = 0, whereas the signs are opposite for all r > 0. 



6. Impurity spectral function 



The impurity spectral function is 
— ImS(w)/7r 



A(w) 



[ui-e d - ReS(w)] 2 + [ImS^)] 5 



(A7) 



If ed 0, then for all r > the spectral function is 
featureless in the vicinity of the impurity energy, 



constant, (A8) 



[ReE( e<i )] 2 + [T(e d )]2 
and takes its low- frequency behavior from T(u>), 

r( w ) 



|w| < |e d |. 



(A9) 



For ed = 0, by contrast, the spectral function exhibits 
nontrivial structure near u> = 0: 



A(u) 



1 



r 2 

1 



7rr r 2 + c;2 ' 



- . , cos — 
7rr(w) 2 



1 



7rr(o;) 
, 7r(gw) 2 ' 



r 7r 



Here. 



.9 = 1 



2r 



(r-l)7rD 



r = 0, 
< r < 1, 

r = l, 

r > 1. 

(A10) 

(All) 



For a flat mixing rate (r = 0), the spectral function con- 
sists of the standard Lorentzian resonance centered on 
oj = 0. For all r > this feature is replaced by a power- 
law cusp, A(uj) ~ Iw/Dl' 1 "''" 1 , such that the spectral 
function diverges for < r < 2 but instead vanishes for 
r > 2. 



c. Conduction-band phase shifts 

The mixing term in the Hamiltonian effectively adds 
one extra state to the band, centered on an energy q 
which is a root of the equation — — ReE(e<j) = 0. 
As a result, each of the original band states is shifted in 
energy from e to e — 5o(e)/ir. Here 5q, the s-wave phase 
shift, satisfies 



at an 



ImE(e) 



e - e d - ReE(e) 



(A12) 



Since one expects the band states to be pushed away from 
the inserted level, the sign of 6o(e) should be opposite to 
that of e — e d — ReE(e). Coupling to the impurity brings 
about a change in the density of states, 



Pimp(e) = S{e- e d ) +ir 1 dS /de. 



(A13) 



For the case of a pure power-law scattering rate with 
(<i =h 0, So(e) w T(e)/e d at low frequencies. The case 
€d = is again more interesting: 



sgn(-e) 



2 r(e 



■ cos — , < r < 1, 



2ln\D/e\ 



9\e\ 



r = l, 



r > 1. 



(A14) 
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Equations (A14) indicate that at e = 0, So jumps through 
(1— rx)n while p; mp has a delta- function peak of weight n. 
In the case r = 0, the standard interpretation of the 
smooth variation of the density of states is that the im- 
purity level becomes completely absorbed into the band. 
This absorption appears to be incomplete for all r > 0, 
with the delta function in pi mp representing the fraction 
of the impurity degree of freedom that remains local- 
ized. For r > 1, not only does the delta function contain 
the entire weight of the original impurity, but a counter- 
intuitive situation arises in which the impurity level re- 
pels delocalized states close to its renormalized position 
(here, Ed = 0) less strongly than it repels states that lie 
further away in energy. 



APPENDIX B: DETAILS OF THE EXPANSION 
OF Fimp AND ximp 

This appendix fills in some of the steps in the deriva- 
tion of the algebraic expressions for impurity thermody- 
namic properties presented in Section [v| In particular, 
we focus on the methods for performing sums over single- 
particle eigenstates arising in the perturbative treatment 
of the discretized effective Hamiltonians introduced in 



Section IV and on extrapolation of the resulting expres- 



sions to the continuum limit. 

Consider expansion of the properties in the vicinity 
of the weak-coupling fixed point. (Analogous arguments 
apply at strong coupling.) The summands encountered 
in these calculations can generally be separated into the 
product of two parts: the first increasing with the index 
j which labels the single-particle eigenvalues, but doing 
so no faster than aL; the second decreasing for large j 
at least as fast as exp(— (3NVj)- Provided that fcgT <C D 
and (3m <§; 1, this decomposition ensures that (a) the 
summand takes its largest value for 1 <C j ^ N, in which 
range the asymptotic forms given in Eqs. (3.7), (3.1C) 
and (3.11) are essentially exact; and (b) the summand 
is sufficiently small for j = 0(1) and for j — 0(N/2) 
that the range of j can safely be extended to run from 
—oo to +oo. If, in addition, A is sufficiently close to unity, 
the sum over j can be well-approximatedEj by an integral 
over the variable u = fint* A 3_l/JV . This procedure, which 
amounts to the replacements 



(JV+l)/2 



i du 

2— 1 In A L a. 



u/Pn, 



f3 N t* 



(2n+l+r)/2 



(Bla) 
(Bib) 
(Blc) 



converts a sum over j to a A-independent integral multi- 
plied by a simple A-dependent prefactor. 

The continuum limit is reached by simultaneously tak- 
ing A — > 1 an d N — > oo in such a manner that (3n [de- 
fined by Eq. (|5.7|)] approaches some value ^« 1. For 



all values r ^ 1, the precise value of (3 drops out of the 
final expression for each thermodynamic property, so this 
prescription produces an unambiguous result. It will be 
shown below, however, that for linear scattering rates the 
leading corrections at the weak-coupling and asymmet- 
ric strong-coupling fixed points depend explicitly on In (3. 
Since (3 has no physical meaning for A — ► 1, the contin- 
uum limit of the discretized thermodynamic calculation 
contains a degree of ambiguity in this special case. 



1. Local-moment regime 



The starting point fo r com puting F- lmp in the local- 



moment regime is Eq. (5.17), which contains four sep- 
arate summations over single-particle eigenstates. Con- 
sider first the unconstrained double sum over j and k. 
One can show that this term represents the second-order 
shift in the ground-state energy of the system due to 
the perturbations Oy and Oj [Eqs. ( |4.7|) ]. This shift is a 
temperature-independent quantity which should not con- 
tribute to the impurity s pecifi c heat. Moreover, since the 
RG transformation [Eq. (2.34)] subtracts off the ground- 
state energy at each iteration, such a term will not be 
detected numerically and can safely be neglected. 

Each of the remaining summations entering -F! mp con- 
tains at least one facto r of pj , which permits appli catio n 
of the transformation (Bl). The last term in Eq. ( [5.17 ), 
which contains a summation over indices j and k ^ j re- 
quires special attention. We find it convenient to define 



^ ) ( t *)2+2r /3 l+, + n 



N 



,,2 ^,2 



A (n-r)JV/2y-' "0j 



a 0j a 0kVjPj 



(Bi) InA 1+r 



— (t*) 1+r (3 r [ ^- 1 A (ri - r '> N ^ 



oo 1+r (N+l)/2 2 

U e« + 1 ^ ((3nvI) 2 -u 2 ' 



(B2) 



For < r < 1, Ep will be dominated by contributions 
from u and k such that ^NVk ~ u — 111 tms case i 

the sum over k can be converted to an integral over v = 
(3 N t*A k -" N , yielding 



Ep = / du dv 



v 



/o Jo v 2 — u 2 e u + l 
Making the change of variables v — > uy, one obtains 

,2r 



(B3) 



E F 



du 
o e 



it - 



dy 



V 



i Jo " v 2 - 1 



for r < 1. (B4) 



The u integral is related to the Riemann zeta function,@ 
while the y integral was evaluated in Ref. |2|. As a result, 
one can write 
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(B5) 



where <fi(x) and tp(r) are denned in Eqs. Q and ( [5~22"1) 
respectively. 



For r > 1, the fc-sum in Eq. ( |B2|) is dominated by the 
largest values of k, and is not well-approximated by an 
integral. If one neglects u 2 in the denominator, the sum 
can be performed directly to give Eq. (B5) once again, 
but with 



ijj(r) 



lnA ( f *Y-i 
1_ A-ir-i)^ 1 i ' 



-(TV + 1) In A, 



r > 1, 
r = 1. 



(B6) 



These manip ulati ons, when combined with Eq. (5.7) 
transform Eq. ( |5.17 ) to 



fim P , „ „r a ai<j)(l + r) ( k B T 

— = 111 z — Oil 



In A (t 

4(^ + ^ 2 ) (gr 



*\2+r 



aD 
1 



l+r 



In A/ (t*) 2 + 2r 
k B T^ L> ' 



In Arc/>(2r) 



aD 



■2%l){r)(j>{rx+r) 



k B T 
aD 



ri+r 



(B7) 



Similar methods can be applied to Eq. Q5.18 ) for Ximp- 
Again, evaluation of the double sum requires the most 
care. We define 



In A 



j x — \ ) (*T +Jr £ 



*\2+2r ol+r+n 



N 



" uV^-1) (N ^ / 2 (p N a 0k ) 2 
o " (^ + 1) 3 £j («) 2 ^ 2 ' 



(B8) 



For r < 1, the fc-sum in Eq. (B8) can be converted to an 
integral, yielding 



f°° v r u 1+r e"(e u -l) 
x ' J v 2 -u 2 (e u + l) 3 



Letting v —* uy, one obtains 

_ fll + r) 
x l + r 



(B9) 



(BIO) 



where </>(x) and are defined in Eqs. (5.5) and ( 5.22 ), 
respectively. For r > 1, direct summation neglecting u 2 



in the denomin ato r gives Eq. (B1C) with ip(r) instead 
defined by Eq. ® . 



The remaining sums in Eq. ( 5. 18] ) are straightforward 
to perform. The resulting expression for the impurity 
susceptibility is 



(.9Ms ) 



J a\ 



(l+r) 



knT 



2 In A (l+r)(t*) 1+r \ aD 



2h a ° ai 



In A (t*) 



(l + r) (k B T 



aD 



*\2+r 



l + r 



4C/ 



(l+r) 



lnA(l + r)(t*) 1 +'' 



k R T 



1 



In Ay (**) 



*\2+2r 



aD 
lnAr0(2r) 

ri+r 



l+2r 



k B T 
aD 



(Bll) 



2r 



+ 2^(r)#r 1 +r) (^\ 



The final step is to extrapolate the expressions for 
i^i mp and Xi mp t o the continuum limit. For In A <C 1, 



Eqs. (2.31), (O), and (3.11) reduce to 



,t*» l + 0(lnA), a 2 n re ±(2n + l+r)lnA. (B12) 



Substituting these values into Eqs. (prf), (B7) and ( Bll|) , 
and then letting In A — * 0, one obtains Eqs. ( 5. 20] ) 
and J5.21D with 



lim ib(r) = (r — 1) 

A->1 



r > 1. 



(B13) 



In t he s pecial case r = 1, application of Eqs. (2.31) 
and (5/7) leads to the result 



lim A/TnA = -21n 

JV-»oo.A^l 



/3fc B T 



(B14) 



where, as stated above, (3 is the limiting value of (3^. 
This in turn leads to the replacement in Eq. (5.26). 



2. Symmetric strong-coupling regime 



The sums entering Eqs. (5.27) and ( 5. 28] ) can also be 
transformed into integrals using the methods described 
in the previous subsection. The leading deviations from 
the fixed-point free-energy and susceptibility arise from 
first-order terms in perturbation theory, so there are no 
double summations to contribute logarithmic corrections 
to the simple power laws in temperature. For small but 
finite In A, one obtains 



k B T 



k B T 



In A (t*A- r / 2 ) 2 - r V aD 



1— r 



(B15) 



and 
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kBTx 



) 2 
-4Ui 



...... = n_ _ 2t~^ 1 ^ 2 ^( 1_r ) 



In A 



In A (t*A-''/2)2-r 

-i 2 



1- 



(l-r)(i*A''/ 2 ) 1 -'- 



aD 

k B T 
aD 



l-r 



(B16) 



where <fi( x ) an d 0(^0 are defined in Eqs. (|5.3|) . Extrapo- 
lation to the conti nuum limit yield s the final expressions 
contained in Eqs. ( |5.30 ) and ( 5.31 ). 
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